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We present a quenclied lattice calculation of the weak nucleon form factors: vector (i<V(f?^)), 
induced tensor {FT{q^)), axial- vector {FA{q^)) and induced pseudo-scalar {Fp{q^)) form factors. 
Our simulations are performed on three different lattice sizes x T — 24^ x 32, 16^ x 32 and 
12'' X 32 with a lattice cutoff of ~ 1.3 GeV and light quark masses down to about 1/4 the 
strange quark mass (m^ « 390 MeV) using a combination of the DBW2 gauge action and domain 
wall fermions. The physical volume of our largest lattice is about (3.6 fm)'', where the finite volume 
effects on form factors become negligible and the lower momentum transfers (g^ « 0.1 GeV^) are 
accessible. The g^-dependences of form factors in the low region are examined. It is found 
that the vector, induced tensor, axial-vector form factors are well described by the dipole form, 
while the induced pseudo-scalar form factor is consistent with pion-pole dominance. We obtain the 
ratio of axial to vector coupling QA/gv = -Fa (0)/-FV (0) = 1.219(38) and the pseudo-scalar coupling 
Fp{0.88ml) = 8.15(54), where the errors are statistical errors only. These values agree 
with experimental values from neutron /3 decay and muon capture on the proton. However, the 
root mean-squared radii of the vector, induced tensor and axial-vector underestimate the known 
experimental values by about 20 %. We also calculate the pseudo-scalar nucleon matrix element in 
order to verify the axial Ward-Takahashi identity in terms of the nucleon matrix elements, which 
may be called as the generalized Goldberger-Treiman relation. 

PACS numbers: ll.15.Ha, 12.38.-t 12.38.Gc 



I. INTRODUCTION 

A comprehensive understanding of hadron structure, especially nucleon structure, based on quantum chromody- 
namics (QCD) is one of our ultimate goals in lattice QCD calculations. The latest lattice calculations of nucleon 
structure have been greatly developed with increasing accuracy So far, large efforts by lattice QCD simulations 
have been mostly devoted to studies of electro-magnetic structure of the nucleon and either unpolarized or polarized 
parton distributions in deep inelastic scattering J[l|, [3, S B IB- However, there are only a few lattice studies to be 
completed for the weak nucleon form factors [6|, |7|, which are associated with weak probes of nucleon structure. In 
this paper, we present results from our intensive study of the nucleon matrix elements of the weak current in quenched 
lattice QCD calculations with domain wall fermions (DWFs). 

Experimentally, weak processes meditated by the weak charged current like neutron beta decay n —^ p -\- e'' -\- De, 
muon capture on the proton /i~ -\-p —^ i'^-\-n or quasi-elastic neutrino scattering D^^ -\-p — > /.t"*" -\-n are mainly exploited 
for studying the weak nucleon form factors, while available information obtained from the experiment of the neutral 
current weak process such as semileptonic elastic scattering u -\-p v ^ p is still limited. The weak current is known 
to be described by a linear combination of the vector and axial-vector currents. In general, four form factors appear in 
the nucleon matrix elements of the weak current. Here, for example, we consider the matrix element for neutron beta 
decay. In this case, the vector and axial- vector currents are given by V^(x) = u{x)^ad{x) and A'^{x) = u{x)^albd{x) 
and then the matrix element is expressed by 

(p|F+(x) + A+{x)\n) = Up [icFviq^) + cj^pqpFTiq^) 

+lc.75FA{q^) + tqc.l5Fp{q^)] u„e*«■^ (1) 

where q — Pn — Pp is the momentum transfer between the proton (p) and neutron (n). The vector (Fy) and induced 
tensor (Ft) form factors are introduced for the vector matrix element, and also the axial-vector (Fa) and induced 
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pseudo-scalar (Fp) form factors for the axial- vector matrix element. The vector part of weak processes are related 
to the nucleon's electro-magnetic form factors, which are well measured up to large momentum transfer by electron 
scattering [1], through an isospin rotation. Based on the conserved- vector-current hypothesis, the vector and induced 
tensor form factors are well understood by knowledge of the electro-magnetic structure of the nucleon. 

In the axial- vector part of the weak process, the axial- vector coupling gA = FA^q^ = 0) is most accurately measured 
by neutron beta decay, where the extremely small momentum transfer is accessible due to a very small mass difference 
of the neutron and proton. The g^-dependence of FA(g^) can be determined by other processes such as quasi-elastic 
neutrino scattering experiments and charged pion electroproduction experiments. It has been observed that the dipole 
form is a good description for low and moderate momentum transfer, <1 GeV^ On the other hand, the induced 
pseudo-scalar form factor Fp{q^) is rather less known experimentally. The main source of information on Fp{q^) stems 
from muon capture. The induced pseudo-scalar coupling, gp — mf^Fp{q'^) evaluated at q^ — 0.88m^, where is the 
muon mass, is measured by ordinary muon capture (OMC) or radiative muon capture (RMC). Although there is some 
discrepancy between the OMC result and the RMC result [1, [l^ , the new precise OMC measurement by the MuCap 
collaboration, which is nearly independent of /i-molecular effect, yields gp = 7.3 ± 1.1 [ll|. Only a few of the other 
q^ data points on Fpijf') are measured in the low region by a single experiment of pion electroproduction [l3 |. 

Theoretically, in the axial part of such weak processes at low energies, one may consider that spontaneous chiral 
symmetry breaking, which is induced by the strong interaction, plays an essential role. In other words, the axial 
structure of the nucleon would be highly connected with the physics of chiral symmetry and its spontaneous breaking, 
which ensures the presence of pseudo Nambu-Goldstone particles such as the pion. This is empirically known as 
the partially conserved axial- vector current (PCAC) hypothesis fs^], where the divergence of the axial- vector current 
is proportional to the pion field. Applying this idea to the axial- vector part of Eq.(IT]), there appears a specific 
relation between the residue of the pion-pole structure in Fp[q^) and the axial- vector coupling gA known as the 
Goldberger-Treiman relation [Tsl . 

There was the long standing disagreement between experiment and lattice calculations about the axial-vector 
couphng gA- However, the RBC Collaboration finally resolved this puzzle using quenched DWF simulations [fllTH- 
DWFs are expected to provide an implementation of lattice fermions with exact chiral symmetry [l^ [13, [3 ■ In the 
limit where the fifth-dimensional extent Lg is taken to infinity, DWFs preserve the axial Ward-Takahashi identity, 
even at a finite lattice spacing [18]. Although not sufficiently large Ls loses the virtues of DWFs, the explicit chiral 
symmetry breaking with moderate sizes of Lg can be attributed to a single universal "residual mass" parameter TOi-cs, 
acting as an additive quark mass in the axial Ward-Takahashi identity as daA'^ « 2(to/ -I- mres)P°' [H, A very 
small value of TOics, which is typically smaller than 10% of the quark mass, is always achieved at a given Lg around 10- 
20 with the help of improved gauge action s 12011 . This fact greatly simplifies the nonperturbative determination of the 
renormalization of quark bilinear currents [2l| . For a calculation of the axial- vector coupling gA , the chiral symmetry 
is very useful because the renormalization factors of local vector and local axial-vector current operators are equal, 
Zy = Za [21| . This means that the ratio of the nucleon axial- vector and vector couplings, gA/gv, calculated on the 
lattice is not renormalized p^ . Therefore, in DWF simulations, the ambiguity in the renormalization of quark currents, 
which is present in other fermions such as Wilson-type fermions, is eliminated. In Ref. J^], gA — 1.212(27) in the 
chiral limit is obtained from quenched DWF simulations. It underestimates the experimental value of 1.2695(29) [S^l 
by less than 5%. It has also shown that there is a significant finite volume effect between the axial-vector couplings 
calculated on lattices with (1.2 fm)"^ and (2.4 fm)'^ volumes. This observation strongly indicates that the axial-vector 
coupling is particularly sensitive to finite volume effects. Subsequently, the LHPC Collaboration has evaluated the 
axial-vector coupling using domain wall valence fermions with improved staggered sea quark configurations with 
physical volume as large as (3.5 fm)-^ and obtained gA = 1.226(84) at the physical pion mass [1^. Its value again 
agrees with experiment within 5%. 

In this paper, we naturally extend the quenched DWF calculation for exploring the axial structure of the nucleon, 
namely the axial- vector form factor and the induced pseudo-scalar form factor at low q^ as well as the electro-magnetic 
structure of the nucleon. Especially, to evaluate the induced pseudo-scalar coupling gp is one of our main targets, since 
no intensive study has been done to determine this particular quantity in lattice QCD. Recall that the induced pseudo- 
scalar form factor is assumed to be dominated by a pion pole, which give rises to very rapid g^-dependence at low 
q^. The larger physical volume, where the lower momentum transfers are accessible, is required. We therefore utilize 
(3.6 fm)'^ volume where the smallest momentum squared (q^ « 0.1 GeV^) is smaller than measured pion mass squared 
(ml > 0.15 GeV^). We also re-examine the finite volume effect on the axial-vector coupling using three different 
volumes, which include (3.6 fm)'^ together with smaller ones (1.8 fm)'^ and (2.4 fm)"^. Furthermore, we calculate the 
nucleon matrix element of the pseudo-scalar density {p\uj5d\n) to check the axial Ward-Takahashi identity in terms 
of the nucleon matrix elements, which may be called as the generalized Goldberger-Treiman relation [24j. 

Our paper is organized as follows. In Section II, we first present a brief introduction of the weak nucleon form 
factors and the status of experimental studies. In Section III, details of our Monte Carlo simulations and some basic 
results are given. We also describes the lattice method for calculating the nucleon form factors. Section IV presents 
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our results of the four weak form factors as well as the pseudo-scalar form factor on lattice with (3.6 fm)^ volume. 
Especially, the q^-dependences of all measured form factors at low are discussed with great interest. At the end of 
this section, we discuss the consequence of the axial Ward-Takahashi identity among the axial- vector form factor, the 
induced pseudo-scalar form factor and the pseudo-scalar form factor. In Section V, we discuss the finite volume effects 
on the form factors using results from three different volumes. Meanwhile, we also check whether approximated forms 
of g^-dependence of form factors, which are observed at low g^, are still valid even in the relatively high region, 
up to at least q^ « 1.0 GeV^, apart from consideration of the finite volume effects. In Section VI, we compare our 
results with previous works. Finally, in Section VII, we summarize the present work and discuss future directions. 



II. WEAK NUCLEON FORM FACTORS AND EXPERIMENTAL STATUS 

In general, the nucleon matrix elements of the weak current are given by a linear combination of the vector and 
axial- vector matrix elements. Here, let us introduce the vector and axial- vector currents, which are expressed in terms 
of the isospin doublet of quark fields "ip = (""i d)'^ 

V^{x) = ^(x)-i^e^{x), (2) 
Alix) = ^{xho.J^t'^^Pix), (3) 

where i° are the SU{2) flavor matrices normalized to obey Tr(t°t^) = Sab- Then, the nucleon matrix elements are 
given by 

{N{P')\J:'^{x)\N{P)) {NiP')\V:{x)+A'^^ix)\NiP)) (4) 
= uJP'){0:^iq)+0^iq))t^u,iP)e^'^-^, (5) 

where q = P — P' is the momentum transfer between the initial (P) state and the final state {P') and N represents 
the nucleon isospin doublet as = {p, n)^ . Four form factors arc needed to describe these matrix elements: the weak 
vector and induced tensor (weak magnetism) form factors for the vector current, 

Ol{q) = lo.Fv{q^) + a^pqpFriq') (6) 
and the weak axial-vector and induced pseudo-scalar form factors for the axial-vector current (s^ 

Oiiq) = 7o75i^A(<z') + z<7«75i^p(9'), (7) 

which are here given in the Euclidean metric convention [55|. Thus, q^ denoted in this paper, which stands for 
Euclidean four-momentum squared, corresponds to the time-like momentum squared as qj^j = —q^ < in Minkowski 
space. 

The weak matrix elements are related to the electro-magnetic matrix elements if the strange contribution is ignored 
under the exact iso-spin symmetry. A simple exercise in SU{2) Lie algebra leads to the following relation between 
the vector part of the weak matrix elements of neutron beta decay and the difference of proton and neutron electro- 
magnetic matrix elements [l^ : 

{p\ujad\n) = {p\uj^u - djad\p) = (pbTb) - {n\jr\n), (8) 

where = ^u^aU — ■^d'jad. This relation gives a connection between the weak vector and induced tensor form 
factors and the iso-vector part of electro-magnetic nucleon form factors 

F^{q') = Fviq'), (9) 
F^{q^) = 2MNFT{q^), (10) 

where Fi (F2 ) denotes the iso-vector combination of the Dirac (Pauli) form factors of the proton and neutron, which 
are defined by 

{n{p')\jT{x)\n{p)) =uAp') + ^0/3^^2^(9')) nJP), (11) 

where Mn denotes the nucleon mass, which is defined as the average of neutron and proton masses, and N represents p 
(proton) or n (neutron). Experimental data from elastic electron- nucleon scattering is usually presented in terms of the 
electric GE{q^) and magnetic Ga/((7^) Sachs form factors which are related to the Dirac and Pauli form factors [sl. [26j 

G%{e) = ' ^F^^ (12) 

GUi') = F^{q^) + F^{q^). (13) 
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Their normalization at = are given by the proton (neutron) charge and magnetic moment [22 |: 

Proton: G^(0) = 1, G^(0) = = +2-792847351(28), . . 

Neutron: G|(0) = 0, G^^(O) = -1.91304273(45). ^ > 

Therefore, one finds Fy(0) = G^(0) - G^(0) = 1 and 2Mi^Ft{Q) = G^^(O) - G5J^(0) - 1 = 3.70589. As for the 
q'^-dependence of the form factors, it is experimentally known that the standard dipole parametrization Gd{(1^) = 
A? /{A? + q^) with A = 0.84 GeV (or A^ = 0.71 GeV^) describes well the magnetic form factors of both the proton 
and neutron and also the electric form factor of the proton, at least, in the low region [2^. Here, the current 
interesting issues of the g^-dependence of the electro-magnetic form factors at higher q^ are beyond the scope of this 
paper. Recent reviews on the experimental situation can be found in Ref. [26|. The slopes of the form factors at 
q^ = Q determine mean-squared radii, which can be related to dipole masses as (rf) = 12 /Mf [i = E oi M) in the 
dipole form Gi{q^) = Gi(0)/(1 -I- q^ /Mf). The experimental values of the electric root mean-squared (rms) radius 
for the proton and the magnetic rms radii of the proton and neutron are compiled in Table HI These rms radii are 
all equal within errors and arc in agreement with the empirical dipole parameter A. On the other hand, the slope of 
the neutron electric form factor G^(g^) is determined with high precision from double-polarization measurements of 
neutron knock-out from a polarized ■^H or ^He target, while only a small deviation from zero is observed for G^(q^) 
at low q^ 26]. Combined with all four of the electric charge and magnetization radii of the proton and neutron, we 
finally evaluate the rms radii for the weak vector form factor and induced tensor form factor as \/\(n7p) = 0.797(4) 
fm and s/ {{rrY) = 0.879(18) fm, which correspond to the dipole masses. My = 0.857(8) GeV and Mr = 0.778(23) 
GeV. See Appendix A for details. 

The axial- vector form factor at zero momentum transfer, namely the axial- vector coupling = Fa{Q), is precisely 
determined by measurements of the beta asymmetry in neutron decay. The value of gA = 1.2695(29) is quoted in the 
2006 PDG [24I. Nevertheless, kinematics of neutron beta decay are quite limited due to a very small mass difference of 
the proton and neutron. Other experimental methods are utilized for determination of the g^-dependence of FA{q^)- 
For this purpose, there are basically two types of experiment, namely quasi-elastic neutrino scattering and charged 
pion electroproduction experiments. The former suffers from severe experimental uncertainties concerning the incident 
neutrino flux and the background subtraction of elastic events, while model-dependent analysis is somewhat inevitable 
for the latter [§]. Both methods reported that the dipole form FA{q^) — Fa(0)/(1 -I- q^ /M\) is a good description 
for low and moderate momentum transfer, q^ < 1 GeV^. The resulting world average of the dipole mass parameter 
Ma is quoted as Ma = 1.026(21) GeV from neutrino scattering or Ma = 1.069(16) GeV from pion electroproduction 
in Ref. As for a small discrepancy between two averages, it has been argued that within heavy-baryon chiral 
perturbation theory (HBChPT) the finite pion mass correction of —0.055 GeV to the latter value may resolve this 
discrepancy Therefore, one can translate the axial dipole mass into the axial rms radius of \/ ((r^)^) = 0.67(1) fm, 
which is consistently obtained from quasi-elastic neutrino scattering experiments and charged pion electroproduction 
experiments 

On the other hand, the induced pseudo-scalar form factor Fp{q^) is less well-known experimentally [10]. The main 
source of information on Fp{q^) stems from ordinary muon capture (OMC) on the proton, /i^ + p ^ n. One 
measures the induced pseudo-scalar coupling gp — mf^Fp{q^) at the specific momentum transfer for the muon capture 
by the proton at rest as q^ = 0.88m^. The induced pseudo-scalar coupling gp is also measured in radiative muon 
capture (RMC), ji^ + p ^ ^ + + n. Before 2006, the Saclay OMC experiment, which was the most recent OMC 
experiment at that time, reported ((7p'^^)saciay, original — 8.7 ± 1.9 [1^. Combining with the older OMC experiments 
including bubble chamber measurements, the world average for OMC is obtained as (gp'^^)oid Ave = 8.79 ± 1.92, 
which is given in Refs. i9] and [2^. Surprisingly, this value is close to the theoretically predicted value by HBChPT, 
gChPT ^ g 26 ±0.16 [Sj]. However, the novel RMC experiment at TRIUMF [H, Ho] is puzzling: their measured 
value of gp'^^ = 12.4 ± 1.0 is quite higher than the theoretical value as is the OMC value as g^^'^ k. lAg'^^'^. 
This disagreement is reduced by reanalysis with the updated /i"*" lifetime [13]. Then, the updated result of the 
Saclay OMC experiment yields (5p'^^)saciavjiDdatcd = 10.6 ± 2.7. Accordingly, the weighted world average for OMC, 
(5'p'^'")updatod Ave. — 10.5±1.8 givcu in Ref. [l^l, is shifted away from the theoretical expected value, while the updated 
average value is in agreement with the RMC result within its error. Indeed, there is a caveat that the ortho-para 
transition rate in /i-molecular Hydrogen, to which either OMC and RMC results are very sensitive, is poorly known 
due to mutually inconsistent results among two experiments [28l Isij and theory [33 ] . Comprehensive reviews of a 
history oi gp have been given in Refs.fg!] and 

Recently, a new OMC experiment has been done by the MuCap Collaboration [ll|. The MuCap result is nearly 
independent of //-molecular effects in contrast with the previous OMC experiments and the RMC experiment. After 
the electro-weak radiative corrections, which were underestimated in the old literature, are correctly taken into 
account [11], the new precise OMC measurement yields 

5^"^"P = 7.3±1.1. (15) 
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Including the new MuCap result and taking into account the electro-weak radiative corrections, the new world average 
of the OMC results becomes (5p'^^)ncwAve = 8.7 ±1.0 [s^. As for other data of Fp{q^), only a few data points 
are measured in the low region by a single experiment of pion electroproduction at threshold [13] ■ These data 
are summarized in Table [ill Three data points from pion electroproduction at threshold are well fitted by the pion- 
pole dominance form, Fp{q^) = 2M]siFA{q^) / {q^ + fn^) [s^], which is also consistent with the value determined by 
the new OMC result at q^ = 0.88m^. Therefore, the pion-pole dominance is confirmed, more or less, through pion 
electroproduction 11^. 

III. SIMULATION DETAILS 

We work in the quenched approximation and use domain wall fermions (DWFs) to compute the nucleon matrix 
elements of the weak current. We generate ensembles of the quenched QCD configuration with the renormalization 
group improved, DBW2 (doubly blocked Wilson in two-dimensional parameter space) gauge action [sl, [s^l at /? = 
&/g^ = 0.87 (a""'^ « 1.3 GeV), where the residual chiral symmetry breaking of domain wall fermions is significantly 
improved with a moderate size of the fifth-dimension such as L.^ = 16 [20| . Indeed, the residual quark mass for 
Lg = 16 is measured as small 10 in lattice units [20|, which is safely negligible compared with the 

input quark masses in our simulations, 0.02 < mf < 0.08. We work with relatively coarse lattice spacing, a ~ 0.15 
fm [5^, which is determined from the p meson mass (20| . 

To study finite volume effects, numerical simulations are performed on three different lattice sizes xT = x 32, 
16^ X 32 and 12^^ x 32. The spatial extents in our study correspond to La ~ 3.6, 2.4 and 1.8 fm. Quark propagators 
are generated for four bare masses rrif = 0.02, 0.04, 0.06 and 0.08 for L = 24 and three bare masses ruf — 0.04, 0.06 
and 0.08 for L — 16 and 12, using DWFs with Lg — 16 and M5 — 1.8. Details of our simulations are summarized in 
Tabic mil In Table IIVI some basic physics results are compiled from Ref . [13] ■ 

The pseudo-scalar meson (pion) masses computed in these calculations are summarized in Table fVl All fitted values 
are obtained from the covariant single cosh fit. It is clear that there is no visible finite- volume effect on the pion 
mass. Measured values for L = 16 are in good agreement with the values found in Ref. [isj . where point-to-box quark 
propagators are used with the mostly same gauge ensembles, while the point-to-gauss-smeared quark propagators are 
utilized in the present study. Our simulated values of the pion mass range from 0.39 GeV to 0.76 GeV. 

A. Nucleon spectra and Dispersion relation 

In order to compute nucleon masses or matrix elements, we define the nucleon (proton) operator as 

xHt,p)^J2^'''"^^-bc[ul{yi,t)CjM2,t)]u,{y3,t) X 0(yi -x)(/)(y2 -x)0(y3 -x), (16) 

X 

where abc and ud have usual meanings as color and fiavor indices. C is the charge conjugation matrix defined as 
C = "ttjy and the superscript T denotes transpose. The superscript S of the nucleon operator x specifies the smearing 
for the quark propagators. In this study, we use two types of source: local source as 4>{'yi — x) = Sij/i — x) and 
Gaussian smeared source. Here we take yi = y2 = = in our calculation. As for the Gaussian smeared source, 
we apply the gauge-invariant Gaussian smearing [4l|, \^ with N = 30, w = 4.35. Details of our choice of smearing 
parameters are described in Ref. [43| . 

We construct two types of the two-point function for the proton. One interpolating operator at the source location 
is constructed from Gaussian smeared quark fields, while the other interpolating operator at the sink location is either 
constructed from local quark fields (denoted LG) or Gaussian smeared ones (denoted GG): 

CsG{t-t,rc,q) - iTr{7'+(x^(t,q)F(isrc,-q)}, (17) 

with S = L or G. The projection operator V+ = can eliminate contributions from the opposite-parity state 

for q^ — [3, [4^. It is rather expensive to make the Gaussian smeared interpolating operator projected onto a 
specific finite momentum at the source location (tgrc)- However, it is sufficient to project only the sink operator onto 
the desired momentum by virtue of momentum conservation. Thus, the quark fields at the source location are not 
projected onto any specific momentum in this calculation. For the momentum at the sink location (tsink), we take all 
possible permutations of the three- momentum q including both positive and negative directions. 

Nucleon masses and energies are computed by using the LG correlators with the five lowest momenta: (0,0,0), 
(1,0,0), (1,1,0), (1,1,1) and (2,0,0) in units of 27r/L. All fitted values, which are obtained from the conventional 
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single exponential fit, for each volume are summarized in Table [V] Next, we examine the dispersion relation of the 
nucleon state in our simulations. The purpose of this examination is two fold: 1) Our analysis should be restricted 
to the lower momenta that do not suffer from large C(a^) errors. 2) The evaluation of the squared four-momentum 
transfer requires precise knowledge of the nucleon energies with finite momentum. The later point can be achieved 
by an estimation of the energy i?(p) with the help of the dispersion relation and the measured nucleon rest mass M^r 
that can be most precisely measured. 

For ruf = 0.04 the measured values of the nucleon energy, which are obtained from L = 24 (open circles), L = 16 
(open squares) and L — 12 (open diamonds), are compared with the relativistic dispersion relation 

Eip) = ^V^ + Ml (18) 

in Fig. [1] with either the naive discrete (continuum-like) momentum pi = ^rii or the lattice discrete momentum 
Pi = sin[^ni] {rii = 0, 1,2, • • •, (L — 1)) for p = {px,py,pz). We observe that the measured energies -^(p) are 
consistent with the estimated values from the relativistic dispersion formula for continuum-like momenta (dashed- 
dotted curve) and lattice momenta (dashed curve) in the range of our admitted momentum except for the largest 
momentum on the lattice with L = 12 as shown in Fig. [T] . The difference between either choice of the discrete 
momentum is mostly comparable to the statistical errors, while differences increase at the larger momentum. To 
restrict ourselves to low region {q-^ < 1 GeV^), we do not use the two largest momenta on the lattice with L = 12 
for the proceeding analysis. Therefore, it is not a concern to choose what type of the discrete momentum in the 
dispersion relation in our current calculation. We simply choose the continuum-like momentum throughout this paper 
and then evaluate the values of the squared four-momentum transfer q^ with the measured rest mass Mn and the 
continuum dispersion relation p^. 



B. Three point correlation functions 

As discussed in the previous section, under the exact iso-spin symmetry (m„ ~ rrid), the SU{2) current algebra 
leads to the following relations [1, [lB| 

{p\V+\n) - 2{p\V^\p), (19) 

{p\At\n) = 2{p\Al\p), (20) 

where = ■^{u-faU — djad) and = ■^{u^a75U — dja'j^d). Thus, we may calculate the weak transition matrix 
elements by the iso- vector proton matrix elements. 

First of all, we define the finite-momentum three-point functions for the relevant components of either the local 
vector current {J^^i^) = u{x)'-fau{x) — d{x)jad{x)) or the local axial-vector current (J7^(a;) — u{x)jal5u{x) — 
d{x)"fa^5d{x)) with the proton interpolating operator x: 

(X(i',p')j'i^(i,q)x(0,-P)> = €{P,P') X f{t,t',E{p),E{p')) + ■■; (21) 

where the initial and final proton states carry fixed momenta p and p' respectively and then the current operator has 
a three-dimensional momentum transfer q = p — p'. Here, Dirac indices have been suppressed. The ellipsis denotes 
excited state contributions which can be ignored in the case of t' — i ^ 1 and t ^ 1. We separate the correlation 
function into two parts: GaiPiP') which is defined as 

Glip^p') = (-Z7 • P' + MN)Oi{q){-ij ■ p + Mjv), (22) 

where 0^{q) corresponds to either Eq.® or Eq.([7|), and the factor f{t, t', E{p), E{p')) which collects all the kinemat- 
ical factors, normalization of states, and time dependence of the correlation function. The trace of GaiPjP') with some 
appropriate projection operator V for specific combinations of F and a yields some linear combination of form factors 
in each F channel. On the other hand, all time dependences of the factor f{t,t',E{p),E{p')) can be eliminated by 
the appropriate ratio of three- and two-point functions 0, 0] 



R(0^ Cf„,«.p',p) 



C'^^(tsink ^src 1 P ) 

where 



C^<=(t.i„k - t, p)C'=G(< - isrc, p')C^^(i«ink - tsrc, p') 



(23) 



C^<=(tsink - t, p')CGG(t _ t,„, p)CLG(t,u.i. - tsrc, p) 

Clait^l) = jTr{7'(x'^(tsi„k,p')J'J(^,q)F(isrc,-p))} , (24) 
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which are calculated by the sequential source method described in Ref . [T5I| : 

In this study, we consider only the case at the rest frame of the final state (p' = 0), which leads to q = p. Therefore, 
the squared four-momentum transfer is given by = 2Mj^{E{q) — Mm)- Nucleon energy E{q) is simply abbreviated 
as hereafter. In this kinematics, g^{p,p') is represented by a simple notation as GaiQ)- Then, the ratio ((23|) gives 
the asymptotic form as a function of the current-operator insertion time t, 



nit) ^ -TriVGliq)} X 



^2M-^E{E + M 



(25) 



in the limit when the Euclidean time separation between all operators is large, isink S> t S> isrc with fixed igrc and 

^sink- 

We choose particular combinations of the projection operator V and the current operator J'^ (T — V or A). We 
consider two types of the projection operator, — V+Jt and V§ — V+j^jz in this study. The latter projection 
operator implies that the z-direction is chosen as the polarized direction. We then obtain some linear combination of 
desired form factors from the projected correlation functions, 



Tr {V'g^iq)} = Mn{E + Mm) [Fv{q^) ~ {E - M^)f^T(g')] 



T^inSYil)} - -ie,j.qjMM[Fviq^) + 2MMFT{q^)], 



(26) 
(27) 



for the vector currents and {i = x,y, z). Similarly, we get 



-Tr{VIQH<l)] 



Mm{E + M^ 



N 



FA{q^)k. 



E + M 



N 



-Fp{q^ 



(28) 



for the axial- vector current jf^ [i = x,y, z). In this calculation, we use at most the four non-zero three-momentum 
transfer q = (n^ = 1, 2, 3, 4). All possible permutations of the lattice momentum including both positive 
and negative directions are properly taken into account. All three-point correlation functions are calculated with a 
source-sink separation of 10 in lattice units, which is the same in the previous DWF calculation of the axial-vector 
coupling gA [l^. For L = 24, we calculate three-point correlation functions with three different sequential sources 
generated with source-sink locations, [isic , ^sink] = [0,10], [10,20], and [20,30] on a given gauge configuration [stI ]. 
while only a single sequential source with [tsrc^sink] — [0, 10] is utilized for L = 12 and L = 16 calculations. 
In Fig. O we plot the dimensionless projected correlators 



A 



lTT{v'gr{q)} 



A)r = -- 



Mm{E + Mn) ' 
2 I iqyMM 



\TT{vig^{q)y 



iqxM^ 



N 



(29) 
(30) 



as a function of the current insertion time slice t ioi nif = 0.04 on the largest volume (L = 24) as typical examples. 
Good plateaus for all squared three-momentum transfer are observed in the middle region between the source and 
sink points. The quoted errors are estimated by a single elimination jack-knife method. The lines plotted in each 
figure represent the average value (solid lines) and their one standard deviations (dashed lines) in the time-slice range 
3<t<7. 

Similarly, Fig. [3] shows and A^ for the axial vector current, which are defined by 



^ iTr{P|t;f(g)} 



Mn{E + Mn) 



1 (h-^^ing^m , YT^{vig^{q)V 



qzqx 



izqy 



(31) 
(32) 



It is worth noting that in the axial-vector channel the z-direction is chosen as the polarized direction in this study. 
Therefore, the longitudinal momentum (q^) dependence explicitly appears in Eq. (|28p . This fact provides two kine- 
matical constraints on determination of the three-point functions in our calculation. First, there are two types of 
kinematics, 7^ and 92 = in the longitudinal component [i — z) of Eq. except for the case of = 3 where 
qz is always non-zero. Secondly, the transverse components {i — x or y) of Eq. (j28p are prevented from vanishing 
by the kinematics only if = 2 and 3, where two components of the momentum including the polarized direction 
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(z-direction) are non-zero. These are the reasons why h-^{qz — 0) and t^^[qz 7^ 0) are separately plotted in Fig. ^ 
and results of h^f^{qz = 0) at = 3 and A;^ at = 1, 4 are missing there. 

Finally, we recall that the lattice local operators J^i.^) {T or A), which are represented as the quark bihnear 
currents, receive finite renormalizations relative to their continuum counterparts. Thus the renormalized form factors 
require some independent estimation of Zr, the renormalization of the quark bilinear currents, [JT^]™'^ = .^ri^/a ]''^**"^''. 
As mentioned previously, good chiral properties of DWFs ensure that the lattice renormalizations of the local currents 
are equal, Zy = Za, up to terms of order 0{a?') in the chiral limit and neglecting explicit chiral symmetry breaking 
due to the moderate size of the fifth-dimensional extent Ls [Hi ■ In this paper, we evaluate Zy at each quark mass from 
the inverse of JV(0) that should be unity in the continuum under the exact SU{2) iso-spin symmetry and multiply 
four weak form factors by this renormalization factor Zy to estimate the renormalized form factors in the chiral limit. 



C. Vector coupling gv and axial- vector coupling qa 

At zero three-momentum transfer |q| = 0, only A^ and are calculable. Then, these directly yield the values 
of {gvY^^^'"" and (3^)^''"''='' respectively. Our resuhs of (.9y)'''"''=^ (ffA)'''"''^'' and their ratio (5^)"" = 
obtained in this calculation are summarized in Table IVIl where old results for L = 8 and L = 16 calculated in 
Ref. dBl are also tabulated. In Fig. m we show the ratios of the axial to the vector coupling {gA/gvY'^^^^'^^ calculated 
for three different volumes as functions of pion mass squared. Clearly, the finite volume effect on {gA/gvY^^^^'^^ can be 
observed. The larger volume results exhibit milder quark mass dependence, while the smallest volume results show a 
slow downward tendency toward the chiral limit away from the experimental point. Therefore, for the largest volume 
results, we simply adopt a linear extrapolation with respect to the pion mass squared to take the chiral limit. We 
obtain the axial-vector coupling = 1.219(38) at the physical point {711-^ = 0.14 GeV). 

We next examine more details of the finite volume effect on the vector and axial-vector couplings separately. 
Combined results from three different spatial sizes, La « 1.8 fm, 2.4 fm and 3.6 fm together with old results from 
La w 1.2 fm, we plot (^y)'^**"^" and (5^)'''**"^° against the spatial lattice size in the physical imit in Fig.[5l The quoted 
errors in the figure represent only the statistical errors, which are obtained by a single elimination jack-knife method. 
The left (right) figure is for lighter (heavier) pion mass. The axial-vector coupling (17^)''^'*"^° shows the significant 
spatial-size dependence, while we do not see any serious finite volume effect on the vector coupling ((7^)''^""^°. This 
indicates that the observed finite volume effect in Fig. |4] stems from that of (g^)''^**"^''. Clearly, it is observed that 
(5yi)''*'*"^° decreases monotonically with decreasing spatial size L. Therefore, we simply utilize the power-law formula 
to estimate the infinite volume limit of the axial-vector coupling as 

^latticc^^^) ^ 5^"''='=(oo) + hL-"" (33) 

with the power three (n = 3). Horizontal lines in figures represent the values in the infinite volume limit and their 
one standard deviations. The values obtained from the largest volume, (3.6 fm)^, are quite close to the values in the 
infinite volume limit. At a glance, the spatial size over 2.5 fm is large enough to accurately calculate the axial-vector 
coupling, at least within the range of our admitted quark mass. 
We finally quote 

gY = 1-219 ± 0.038(stat) ± 0.024(norm) ± 0.002(volume), (34) 

where the second error is evaluated from a 2 % error stemming from Zy 7^ Za, which was observed previously [T^ 
and the third error is estimated from a difference between the extrapolated value in the infinite volume limit and the 
largest volume result at m/ — 0.04. 



IV. RESULTS OF NUCLEON FORM FACTORS ON A (3.6 fm)^ VOLUME 

In this section, we focus on the results obtained from lattice size 24'^ x 32, which corresponds to physical volume, 
V ~ (3.6fm)'^. The lower momentum is admitted by the larger spatial extent L. Therefore, we can make the shorter 
extrapolation with respect to toward the forward limit, q^ = 0, for nucleon form factors, Friq'^) and Fp{q'^), of 
which values at = cannot be accessible directly due to the kinematical constraint as described before. We also 
discuss the finite size effect on the nucleon form factors, which may be sensitive to the nucleon "wave function" or 
the nucleon "size" squeezed due to the finite spatial extent of lattice volume. In the previous section, the spatial 
lattice-size dependence of the axial-vector coupling shows that spatial lattice size La w 3.6 fm is large enough to avoid 
significant finite volume effect on gA- 
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A. Vector channel 



In the case if spatial momentum transfer q is non-zero, all three-point correlation functions defined in Eas. (|26p and 
(P7)) are calculable. Two independent form factors Fv(q^) and F^iq^) are obtained by 

2n 2AfAr y E - Mn y 

Friq^) = ^(-Ao^ + A^) (36) 

at finite q^. 



1. Dirac form factor Fv{q^) 



First, we show quark mass dependence of the Dirac form factor Fylq"^). In Fig. [51 we plot the normalized Fv(g^) 
by -FV(O) as a function of four-momentum squared q^. Different symbols represent the values obtained from different 
quark masses nif. There is no large my-dependence, while it seems that the smaller quark mass makes the q^- 
dependence steeper. The Dirac form factor is supposed to be the dipole form at low q^: 



FyjO) 

(1 + qyM^y 



Fvil') = 7T^6^.^ (37) 



where My denotes the dipole mass. A dashed curve in Fig. [6] corresponds to the dipole form with the empirical value 
of the dipole mass My — 0.857(8) GeV, which is evaluated from the electric charge and magnetization radii of the 
proton and neutron as described in Appendix A. 

In order to see how our measured Fv{q^) has an expected dipole form, we define the following quantity only at 
nonzero momentum: 



which should provide independent plateau if the g^-dependence of Fv{q^) ensures the dipole form. We call this 
quantity the effective dipole mass hereafter. In Fig. [71 we show the effective dipole-mass plot for the Dirac form 
factor at m/ = 0.02 as a typical example. Horizontal solid and dashed lines represent the fitted dipole mass obtained 
from a correlated fit to Fy(g^) using the dipole form (|571) and its one standard deviation. The dotted lines shows 
the empirical dipole mass. My = 0.857(8) GeV. Clearly, there is no appreciable g^-dependence of the effective dipole 
mass within statistical errors. Even at the highest q^ w 0.44 GeV^, the fitted dipole mass agrees with the value of the 
effective dipole mass. Therefore, we conclude that the dipole form describes well the g^-dependence of our measured 
Dirac form factor Fv{q^). This observation is consistent with previous studies [1, |4|. 

Fig. [8] shows the quark-mass dependence of the fitted dipole mass. As seen from this figure, the quark-mass 
dependence is rather mild and then there is no appreciable curvature as a function of the pion mass squared. Therefore, 
we simply adopt a linear extrapolation with respect to the pion mass squared to evaluate the value of the dipole mass 
oi Fv{q'^) in the chiral limit. Diamond symbols in Fig.[8]are extrapolated values for the chiral limit (m^ = 0) and the 
physical point (m,r = 0.14 GeV) and the solid line represents the fitted line. Our measured dipole masses of the Dirac 
form factor are much larger than the experimental value. Here we recall that the root mean-squared (rms) radii can 
be determined with the corresponding dipole mass as yj{r^) = Vvl/Mv- The larger dipole mass My means that the 
spatial size of the nucleon in our simulations is smaller than the physical one. This may be attributed to the missing 
large "pion-cloud" contribution since it is well known that the mean-squared radius {vy) receives a large pion loop 
correction, which leads to a logarithmic divergence in the chiral limit in heavy baryon chiral perturbation theory [46l |. 
Indeed, the present calculation is still far from the chiral regime: our smallest pion mass is around 0.39 GeV which 
is comparable to the lightest pion mass in the most recent lattice study of nucleon electro-magnetic form factors [3] ■ 
Although the expected chiral behavior is not guaranteed in the quenched approximation, in the present study the 
estimation of systematic errors stemming from quenching and a long chiral extrapolation is beyond the scope of this 
paper. Rather we would like to see how large volume can be fitted for studying the structure of the nucleon, namely, 
the nucleon form factors without significant finite volume effect. Studies of the finite size effect on the nucleon form 
factors using results obtained from three different lattice sizes will be presented in the next section. 
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2. Pauli form factor Ft (q^) 

In Fig. [9l we show the Pauh form factor F2{q^) = 2MNFT{q^) as a function of four- momentum squared q^. 
The form factor plotted here is scaled by the renormalization factor Zy = 1/Fy(0) to get the renormalized one, 
^2°"(9^) — ^y-p2('?^)- In contrast with Fig. [SI large my-dependence is observed. As well as the Dirac form factor, the 
Pauli form factor is phenomenologically supposed to be the dipole form at low q^: 



where the value of at q^ ~ is associated with the difference of the proton and neutron magnetic moments, 

lip — Hn = 1 + F2"^{0). This dipole form is commonly adopted as a fitting form of the q^ extrapolation to evaluate 
i^™'^(0) in published works [1, We also plot the dipole form with the empirical value of the Pauli dipole mass 
Mt = 0.778(23) GeV and the experimental values of /Xp and /x„ in the same figure. Our results of gradually 
approach this dipole form as to/ decreases. Indeed, data points for ruf = 0.02 in the range of our calculated q^ 
follow the experimental curve within the statistical error. However, if we apply the dipole form to the data, our 
obtained values of i^|°"(0) from dipole fits are somewhat underestimated in comparison with the experimental value 
— A'n — 1 = 3.70589 as shown in Fig. [TUl 

In contrast to the Dirac form factor, we cannot try the effective dipole mass plot for a justification of the applied 
dipole form, since we do not have data of ^^'^"(O) without the q^ extrapolation. Instead, we consider an independent 
observation for the difference of fip and which can be derived from the forward limit of the ratio of the magnetic 
form factor GM{q'^) and the electric form factor Gsiq'^)- The ratio is calculated by a different combination of and 
A^ [13 as 

GMjq') ^ G^(g^) ^ A£ 

GE{q^) G'r(g2) AJ" ^ ^ 

Experimentally, it is known that this ratio shows nog^-dependence at low q^ since both form factors are well fitted by 
the dipole form with the comparable dipole masses [liHl]. Therefore, this ratio may yield the constant value identified 
to Hp — Hn — G™j"(0). Indeed, in our calculation, the proposed ratio ([^0]) exhibits no appreciable g^-dependence in 
the range of our calculated q^ . We may use a simple linear fitting form with respect to four- momentum squared 
for an alternative evaluation of the value /ip — /i„. Fig. [10] shows that two determinations to evaluate /ip — /i„ — 1 are 
consistent with each other. 

In Fig.[TTl we plot values of ^p — /i„, which are evaluated by two determinations, 1-|-F2°'^(0) and Gm(0)/G_b(0), as 
a function of pion mass squared toJ. As described above, both determinations fairly agree with each other. Although 
the values measured at two heaviest points are consistent with the experimental one, a strong toJ dependence appears 
near the chiral limit and then the extrapolated value tends to somewhat underestimate the experimental data. Here, 
we simply adopt a linear fit with respect to to^ regardless of the fact that a slight downward curvature is observed in 

We also extrapolate the Pauli dipole mass Mt to the chiral limit in Fig. [T^l Again, we use a simple linear fitting form 
for the chiral extrapolation. The value obtained at the physical point is about a 20% overestimation in comparison 
with the experimental one, the same as in the case of the Dirac dipole mass. This indicates that corresponding rms 
radii are somewhat smaller than the actual nucleon size. Finally, all fitted results with the dipole form for both form 
factors and their extrapolated values to the chiral limit are summarized in Table \VTl[ 

B. Axial- vector channel 

In the axial- vector channel, two independent form factors FA{q'^) and Fp{q^) can be evaluated separately by 

FA{q^) = Af (g, = 0), (41) 
Fp{q^) = A^/Mn (42) 

at finite q^. It should be reminded that Af^{qz = 0) at = 3 and A^ at = 1, 4 are not obtained directly from 
corresponding three-point functions due to the kinematics as described in the previous section. However, instead, we 
can evaluate them by using a relation 

Atiq, ^ 0) = Atiq, = 0) - , A^, (43) 

where Af (g^ ^ 0) are always calculable at finite g^. 
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1. Axial-vector form factor FAi'f') 

Fig. 1131 shows quark mass dependence of the axial- vector form factor .Fa ((7^). The vertical axis is normalized by 
Fa[0) and the horizontal axis denotes the four-momentum squared in physical units. Different symbols represent 
the values obtained from different quark mass m f . The axial- vector form factor is phenomenologically fitted with the 
dipole form, at least at low q^, as well as the Dirac and Pauli form factors [9]: 



Fa{0) 



FA{q') = , 1\ (44) 



where Ma denotes the axial dipole mass. A dashed curve in Fig. [13] shows the dipole form with an experimental value 
of the axial dipole mass Ma = 1.026(21) GeV 0. There is a similarity here in comparison with Fig. [6l No large 
m/-dependence is observed. Even at the smallest quark mass ruf — 0.02, where the corresponding pion mass is less 
than 400 MeV, our observed FA{q^) is far from the experimental curve. Indeed, the axial-vector form factor FA{q^) 
is flatter than the experimental one, similar to what we observe in Fv{q^) and ^2(9^) = '^M^Fxiq'^)- This again 
indicates that the nucleon size in coordinate space shrinks away. The similar observation is reported in Ref. 

Next, to see how the dipole form is fitted to our measured Fa((?^), we show the effective dipole mass plot, which 
was defined similarly to Eg. ([55)1 . Fig. [T3]is plotted for m/ = 0.02 as a typical example. We also include the fitted 
Ma, which is obtained from a correlated fit to Fa (5^) using the dipole form (|^^ . with its one standard deviation as 
solid and dashed horizontal lines. All momentum points except the third one, which deviates from the fitted value 
by about 2(t, are located within the horizontal lines. Here, we remark that the third momentum point of Fa (9^) at 
m/ = 0.02 in Fig. [13] is slightly dropped from the values measured at other quark masses. We then stress that the 
case of TO/ = 0.02 is the worst example. Indeed, it is found that the effective dipole mass plot for heavier to/ are 
quite consistent with the fitted values in all q^ range that we measured. From this observation, we conclude that the 
q'^-dependence of our measured FA(g^) can be well described by the dipole form (|44| in the range of our utihzed g^, 
(q^ < 0.44 GeV^). This is quite consistent with the phenomenological knowledge on the g^-dependence of FA(g^). 

We show the quark-mass dependence of the fitted axial dipole mass as a function of the pion mass squared in 
Fig. [TS] All measured values are listed in Table IVIIII We find that the quark-mass dependence is somewhat milder 
than the dipole masses for the Dirac and Pauli form factors. Clearly, there is no appreciable curvature as a function 
of the pion mass squared. As before, we simply adopt a linear extrapolation for the axial dipole mass Ma toward the 
chiral limit. The extrapolated values (diamond symbols) overestimate the experimental one marked by the asterisk 
in Fig. \TE 

As pointed out before, our observed "size" of the nucleon in coordinate space is much smaller than the experimental 
one. We will see that there is no significant finite volume effect, which may cause the "size" of the nucleon to be 
squeezed on the lattice with (3.6 fm)'^ volume. Thus, this observed "squeezing", which is evident from the broadening 
of the form factors, may be attributed to the missing contribution of the "pion-cloud" surrounding the nucleon outside 
of the chiral regime. Interestingly, however, the ratio of the axial dipole mass to the Dirac dipole mass is in very good 
agreement with the experiment. In Fig]!!] we show the ratio Ma/ My together with the ratio Mx/My as a function 
of the pion mass squared. The quark mass dependence of both ratios is found to be very mild in our observed range 
of TO^. All measured values of Ma/ My and Mx/My, which are listed in Table llXl are fairly comparable to their 
respective experimental values. We obtain M^/My = 1.285(73) and Mx/My = 0.869(57) at the physical point by 
using a simple linear extrapolation. 



2. Induced pseudo-scalar form factor Fp{q^) 

First, we show the quark- mass dependence of the induced pseudo-scalar form factor Fp{q^) in Fig. 1171 In contrast to 
the axial-vector form factor Fa((7^), significant TO/-dependence is observed especially in the lower q^ region (g^ < 0.3 
GeV^). This might be associated with the pion-pole contribution to Fp{q^), which is expected theoretically. Indeed, 
the partially conserved axial-vector current (PCAC) hypothesis and pion-pole dominance (PPD) predict that the 
induced pseudo-scalar form factor approximately behaves like 

^PPD(^2) ^ 2Af^.Fr(g^) ^ (45) 
g^ + TO^ 

which becomes exact in the chiral limit where the pion is massless {m^ ~ 0) 0, [s^. The single pion electroproduction 
experiment also supports the PPD form J^]. Here, to see how the pion-pole behavior is preserved in Fp[q^) measured 
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in the quenched calculation, we consider the following ratio 

(46) 



^PPD 



FPPD(92)' 



which is inspired by the above PC AC prediction. If the measured Fp{q^) has exactly the same form described in 
Eq. (|45|) . the ratio app^ yields the value of unity in the entire region. 
In Fig. [THl we plot the above defined ratio 

'-^PPD ^ 

function of four-momentum squared q^ . This figure shows 
two important features. The significant quark-mass dependence observed in Fig. 1171 almost disappears as expected. 
Furthermore, there is no appreciable g^-dependence in ctpp^ . Clearly, four different q^ points of a^^^ reveal q^ 
independent plateau within the statistical errors. We simply take the weighted average of ctppo within four measured 
q^ points, then plot them against the pion mass squared. As shown in Fig. the average values of Ofpp^ gradually 
approach unity as the pion mass decreases. However, a simple linear extrapolation of app^, yields a value slightly 
smaller than 1 even in the chiral limit. As a result, measured F^p^{q^) is quit well described by the PPD form with a 
multiplicative (quenching) factor app^ < 1. 

F|f"(g2)«app„ xFr°(9'). (47) 

The validity of the PPD form is also tested by the other analysis. Following the analysis done in the previous study of 
the g^-dependence of Fp{q^) 0, we apply the monopole fit to the ratio Fp{q^) / F^iq^)- The satisfactory consistency 
between the fitted monopole mass and the measured pion mass is observed in our DWF calculation [5gl. while both 
quenched and unquenched Wilson simulations fail to exhibit the correct pion-pole structure of Fp{q'^) [7|]. 

Next, we evaluate the induced pseudo-scalar couphng, which is defined by (gp)"^" = m^Fp'^{'d.SSm?^) where to^ is 
the muon rest mass and F'^p'^{q^) = Fp{q^)/Fv{Q)- The specific momentum transfer for muon capture {q^ = 0.88to^) 

is still far from our lowest momentum transfer (g^ w 0.1 GeV^) so that the determination of (gp)"^" requires the 
extrapolation of Fp^'^{q'^). We have already learned that the q^-dependence of measured Fp'^{q'^) is well described by 
the PPD-likc form (|47p in the low q^ region. Therefore, the induced pseudo-scalar coupling can be evaluated by 

<"'■>"" ^ ' -„.FT'(OM„i), (48) 

where F^''"(0.88m^) is precisely determined through the q^ interpolation with the help of the dipole form and very 
accurate data F^'^{0). The pion mass in Eq. (j48|) is simply replaced by its physical value in order to subtract the 
dominant source of the large m/-dependence. In Fig. [20l we plot the resulting value of (ffp)"^"" (square symbols) as 
a function of the pion mass squared. All measured values are listed in Table lylTTl Although there still remains the 
explicit dependence of the quark mass, the simple linear extrapolation yields {gpY°'^ — 8.15 ± 0.54 at the physical 
point (m^ = 0.14 GeV). 

Here, the observed TO/-dependence stems from that of measured Mn, since both ttpp^ and FA^{q'^) have a very 
mild quark mass dependence. To diminish the explicit my-dependence, we may evaluate the dimensionless prefactor 
in Eq.dUl) with the experimental values of = 105.7 MeV, = 139.6 MeV and Mn = 938.9 MeV instead of 
using measured values. We then obtain a more simple form as (gp)"" = 6.77 x appp,F2'"(0.88m^), which is similar to 
the known phenomenological form for gp beside the quenching factor ttpp^ [idj . Indeed, evaluated values using this 
simple formula have no appreciable m / dependence as shown in Fig. [20l After the linear extrapolation, we obtained 
(Sp)'^™ = '^•31 ± 0.39 at the physical point. Two determinations provide consistent results within their statistical 
errors. Of course, the latter determination is rather phenomenological, then we prefer to quote the former value for 
our final result. We finally quote 

(5p)™" = 8.15 ± 0.54(stat) ± 0.16(norm), (49) 

where the second error is an estimate of a 2% error stemming from Zy 7^ Za, the same as in the axial- vector coupling 
9A- 

This value is to be compared with the most recent experimental value g'^p^ — 7.3 ± 1.1 from the MuCap experi- 
ment [nij where the obtained value of g°p^ is nearly independent of /i-molecular effects. We also quote the prediction of 
chiral perturbation theory, g'p^^'^ = 8.26±0.23 and the new world average of experimental values, g°p^ = 8.7±1.0 [1^ 
obtained from ordinary muon capture including the new MuCap result. 

Phenomenologically, the residue of the pion pole in the induced pseudo-scalar form factor is related to the pion- 
nucleon coupling g^^N n [35l | . The induced pseudo-scalar form factor should be expressed as 
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near the pion pole [q^ « ~'^^) [g. [Sal with the renormaHzed pion decay constant i^^, which is defined as 
ZA{^\daA'^{x)\'iyb{(l)) ~ m'^F^^Sabe'''^'^ [5§|. This parameterization provides a way to evaluate the pion-nucleon cou- 
pling g-nNN from the measured induced pseudo-scalar form factor as follows: 

g^^^ = ^hm [q^+ml)^^-^ (51) 
X MNFn-ml), (52) 



where the second equality follows from our observed form (|47)) on Fp'^{q^). The value of F^'^ {—ni^) is evaluated by 
the dipole form with measured ^^^'"(O), Ma and TOtt- We then obtain 

g^^jv = 10.4± l.O(stat) (53) 

at the physical point. Our obtained value is about 20% smaller that a recent estimation q-^nn = 13.32 ± 0.09 
((7jjY^/47r — 14.11 ± 0.20) obtained from forward TriV scattering data (irj . 

C. Pseudo-scalar channel 

1. Pseudo-scalar form factor Gp{q^) 

In this study, we also calculate the pseudo-scalar nucleon matrix element 

{N{P')\P^ix)\NiP))=uAP')0^iqruAP)e''-\ (54) 

which is associated with the axial- vector matrix element through the axial Ward-Takahashi identity. Here, P°(x) is 
a local pseudo-scalar density, P°(a;) = 'i[){x)^^t°"ip{x). The pseudo-scalar matrix element can be described only by a 
single form factor, which is called the pseudo-scalar form factor Gp{q^): 

OP{q)^^,Gp{q^). (55) 

To extract the form factor, therefore, we simply calculate the following trace of G^{q), which represents the spinor 
structure of the corresponding three-point function, with the projection operator Pf: 

^Tr {P|e^(g)} = iq,MNGp{q^), (56) 

where the definition of G^{q) is given by Eqs. (I^T]) and with the local current J^{x) — u{x)j5u{x) — d(x)'^^d(x) . 
It is apparent that non-zero three-momentum q 7^ is required to access the pseudo-scalar form factor. In other 
words, Gp{q^) in the vicinity of 5^ = cannot be evaluated without q^ extrapolation. 

In Fig. [211 we show the TO/-dependence of the pseudo-scalar form factor Gp{q^). Significant TO/-dependence is 
observed in the lower q^ region, similar to the induced pseudo-scalar form factor Fp{q^). We will discuss the g^- 
dependence on Gp{q^) from the viewpoint of pion-pole dominance later. 

2. Test for the axial Ward- Takahashi identity 

First, we address the question whether our domain wall fermion (DWF) calculations of nucleon form factors satisfy 
the axial Ward-Takahashi identity. A similar study has been recently done with the Wilson fermions in both quenched 
and unqucnched simulations 

In the limit where the fifth-dimensional extent Ls is taken to infinity, domain wall fermions preserve the axial 
Ward-Takahashi identity at nonzero lattice spacing [l^. The axial Ward-Takahashi identity for the DWF with SU{2) 
iso-spin symmetry is 

d^Alix) = 2mfJ^{x) + 2J»,(x), (57) 

where is the partially-conserved axial- vector current, which is point split and requires sums over the extra fifth 
dimension of the DWF, Jg is a usual bilinear pseudo-scalar density corresponding to P", and Jg^ is a similar pseudo- 
scalar density defined at the midpoint of the fifth dimension. The "midpoint" term Jg^ is responsible for the explicit 
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chiral symmetry breaking due to the finiteness of the fifth-dimension fl§\. With moderate Lg, this effect can be 
described by the so-called residual mass term m^cs- Then, Eq. ((57|) can be approximately represented by 



dc^A^Jx) « 2(m/ + m,.cs)P"(a;). (58) 

This residual mass m^-cs is determined by (0| J5q|7r)/(0| J5 |7r). The value of mj-cs is known to be small in this calculation. 
(See Table lYl) 

For a practical reason, we did not use the conserved axial-vector current for evaluating the nucleon axial matrix 
element in this study. Instead, we use the local axial-vector current = V'TmTs^^V'j which may be related to the 
conserved axial-vector current as Af_i = ZAAf^ + 0{a^, mjo^). It is worth mentioning that Zy = Za is satisfied up to 
small discretization errors of 0{a^) in the chiral limit [15| . We also did not measure the nucleon matrix element of 
J§g, therefore we cannot fully check the axial Ward-Takahashi identity (AWT) in terms of the nucleon matrix element 
in this paper. Instead, we test the following ratio, which may have no apparent q^-dependence. 



(59) 



2m/Gp(g2) 

This ratio ([55)1 is associated with the following identity [60j : 

ZA{N\da^A''Jx)\N) =2m^^^{N\P''ix)\N), (60) 

where rrij^^,^ = aj^-^,j.mf, which is expected to be comparable to mf + irises up to terms of order 0{a^, m/a^). 

As shown in Fig |22i indeed, there is no appreciable g'^-dependence in the ratio a^^^^ for each m/. Four different 
points of q;^„^ reveal a g^-independent plateau within the statistical errors. We evaluate the weighted average of 
<*AWT by using all four measured points. The obtained values of a^wx tabulated in Table 1x1 Deviation from 
unity is getting large as mf decreases [6l| . This indicates that a^wx possess a 1/m/ term, which is induced by 
the presence of the additive mass shift in the axial Ward-Takahashi identity such as mres- To see this point clearly, we 
plot the modified ratio as m/(a^„^ — 1), which can be interpreted as the difference between m^^^ and m/. Fig. [23] 
shows modified ratios as functions of four- momentum squared q^ for all four values oi mf. There is again no visible 
q^ dependence. Moreover, all mf results are consistent with each other within statistical errors, mghift = f^Awi ""^/j 
which is given by taking the weighted average of four q^ points in Fig. I23i corresponds to the relative amplitude of 
{N\J?^^\N) to the usual pseudo-scalar matrix element {N\JS^\N). More precisely, TOshift is expressed by 

w,hift= (jv|Jg|jV) ^ ' 

Therefore, supposing that {N\J^^\N) / {N\J^\N) sa (0| Jj?^|7r)/(0| Jj?|7r), we expect Wshift ~ TOios besides 0{a^,mfa?) 
corrections. We plot mghift against the pion mass squared in Fig. 1241 The mild wy-dependence allows us to take 
a linear extrapolation for mshift to the chiral limit. At mt = 0, we obtain TOshift = 0.0073(12), which is about one 
order of magnitude larger than TOj-os = 5.69(26) x 10"'* [23| contrary to our naive expectation. A few % level 0{a^) 
correction, which is observed in the difference between Zy and Za cannot account for this discrepancy. To resolve it, 
it is necessary to calculate the relative amplitude of {N\J^^\N) to the usual pseudo-scalar matrix element {N\J?^\N) 
directly We plan to study {N\J^^\N) as well as {N\A'^\N) with the conserved axial- vector current Aa in our extended 
work 14811. 



3. Test for the pion-pole dominance on Gp{(f) 

According to the pion-pole dominance of Fp{q'^), we may expect that the pion-pole dominance holds even in Gp{q^). 
As described in Appendix B, a naive pion-pole dominance hypothesis predicts the ratio of the pseudo-scalar form 
factor and induced pseudo-scalar form factor, will not depend on q^ at low q^ but will exhibit a constant value, related 
to the low energy constant Bq. Indeed, this is not the case. In Fig. [25l we show the ratio of our measured Gp(g^) 
and Fp{q^) as a function of momentum squared q^ . There is a linear-like q^ dependence, which strengthens as the 
quark mass decreases. However, as we will describe below, we confirm that the pion-pole dominance on Gp{q^) still 
remains valid in our calculation. 

What we have observed in the previous subsection can be interpreted as a consequence of the axial Ward-Takahashi 
identity among three nucleon form factors: 

2M^Fr (g2) - g2^^|,-(g2) « 2m,^,Gp{q^). (62) 
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Combined with this relation and the important observation of pion-pole dominance in Fp{q'^) fEa. (|T7)l ). one may ex- 
pect that the g^-dependcncc of Gp{q^) is mostly described by the pion-pole dominance form with a slight modification, 
which corresponds to an extra g^-dependence caused by the fact that ctpp^ 1: 

1 + (1 - 

'-^AWT 

where the naive pion-pole dominance form [49l | is defined as 

2m^GPP°(<z2) = 2M^Fr(g') Vt^- (64) 



g2. 



This residual g^-dependence due to app^, 7^ 1 is supposed to be responsible for the g^-dependence in the ratio of our 
measured Gp{q^) and Fp™(g^) 



F-ril') "„„"™, F^^CO^) 2m 



AWT 



where Appd(5'^) = (1 -I- (1 — ot^^^)^) / a^^^. It is clear that the residual g^-dependence becomes large as goes 
to zero. This feature is in agreement with what we observe in Fig |25l If we multiply the ratio by the corresponding 
factor AppD(g^), which is evaluated with measured ctppp, and m^, the linear- like g^-dependence indeed disappears as 
indicated in Fig |26l Four different points of this ratio reveal q^ independent plateau within their statistical errors. 
We then can evaluate the weighted average by using all four measured q^ points to obtain the value of Gp{q^)/F'^p^{<f') 
in the limit of q^ — > 0, which corresponds to m^/(2TO^-^^^appp,). 

Although G'p(0)/Fp°"(0) is associated with the bare value of the low-energy constant in the pion-pole dominance 
model as discussed in Appendix B, this prediction is slightly modified as 

lini = ^J^11±^B,^ (66) 

where the low-energy constant Bq is defined by the relation — 2{mf + mycs)Bo. Because of the fact that 
TTires 7^ "T-shift , the valucs of G p (0) / Fp'^ (0) X Qfppp, deviate from Bq as shown in Fig. [571 The horizontal solid and 
dashed lines represent the reference value of Bq and its standard deviation, which are evaluated from fitting squared 
pion masses to the linear function cq + Ci ■ rrif. The fit yields the low-energy constant Bq = ci/2 = 2.0705(93) in 
lattice units. The ratio of fitting parameters cq and ci gives rise to the value of 1.75(15) x 10^'^, which slightly 
overestimates TOios quoted in Ref. [20|. This is simply because our fitting is performed in the rather heavier quark 
mass region (0.02 < to/ < 0.08). For comparison, values of toJ/2to/ and to^/2(to/ -t-TOj-os) are also plotted as square 
and circle symbols in Fig. 1271 This observation may raise a question whether the single universal "residual mass" 
parameter exists. However, what we observe here is very sensitive to the correct chiral behavior of the nucleon matrix 
elements. The quenched approximation may provide unknown quenching effects in nucleon matrix elements near the 
chiral limit. In this sense, dynamical simulations are much preferable to investigate it further. 



V. FINITE VOLUME EFFECT ON NUCLEON FORM FACTORS 



As we discussed in Sec. IIII Cl we have found a significant finite volume effect on the axial- vector coupling qa, while 
there is no appreciable effect on the vector-coupling gy- In this section, we test for finite volume effects on all of the 
five nucleon form factors computed in this study. Unlike those couplings gv and gA, which are defined at = 0, it 
is hard to compare values of the form factor at non-zero q^ among different spatial sizes L. This is simply because 
non-zero q^ values are discrete in units of (27r/L)^. In this study, however, our largest spatial size (L = 24) is twice as 
large as the smallest one (L = 12). The q^ value for q = ^(2, 0, 0) at i = 24 coincides with the one for q = ^(1, 0, 0) 
a.t L = 12. Therefore, at least, a single value of non-zero q^ is common between two different lattice sizes. 

In Fig. [251 we show results for the vector form factor Fyiq"^) for three spatial sizes. The left (right) panel is for 
TO/ = 0.04 (to/ = 0.08). Different symbols represent the values obtained from simulations with different spatial lattice 
sizes. Solid curves are dipole form fits results on the largest volume (L — 24) as we described in Sec lIV ATI These 
curves should be capable of exposing finite volume effects on the form factor. Data points from smaller lattice sizes of 
either i = 16 or L = 12 at heavier to/ agree well with the dipole fits, while results obtained from the smallest lattice 
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(L = 12) at lighter m/ seem to be slightly away from them. However, the values obtained from L = 24 and L = 12 
at « 0.44 GeV are not significantly different in the statistical sense. 

We next show the same type of figures for the induced tensor form factor Fxiq"^) in Fig.[29l Qualitative features are 
quite similar to the case of the vector form factor. In the case of the heavier quark mass (m/ = 0.08), all data points 
follow the solid curve fairly well, which is the dipole form fitted to the data of the largest volume (L = 24). On the 
other hand, at lighter quark mass {mf = 0.04), the data points obtained from the smallest lattice volume (L = 12) 
slightly overestimate the solid curve. Again, the value for L = 12 at « 0.44 GeV^ is not significantly away from 
the value for L = 24 in the statistical sense. 

Next, let us examine two form factors in the axial-vector channel, where the axial-vector coupling gA significantly 
suffers from the finite volume effect. In Fig. [30l the axial-vector form factor _Fa((7^)/Fa(0) for three spatial sizes is 
shown. Gross features are mostly similar to both the vector form factor and the induced tensor form factor. At the 
heavier quark mass, it is observed that all data points follow the solid curve fitted to data of the largest volume with 
the dipole form Although two values at « 0.44 GeV^ obtained from both largest and smallest volumes agree 
each other within statistical errors, the finite volume effect seems to be non-negligible in the lighter quark mass region. 
On the other hand, the finite volume effects do not show up in the induced pseudo-scalar form factor except for the 
lowest q^ value for L = 16, as indicated in Fig. [311 The solid curves are obtained by fitting data of the largest volume 
with the PPD-like form (|47p. where the pion-pole structure is essential at low q^. At this moment, we do not have 
any explanation why the lowest q^ data for L = 16 deviates from the solid curve. We also show the form factor in the 
pseudo-scalar channel, namely Gp{q^), for three lattice sizes in Fig. [321 At either quark mass nij — 0.04 or 0.08, all 
data points agree well with the solid curves fitted to data of the largest volume by the modified PPD form All 
data that appear in Fig. [28ll32l are tabulated in Tables IXHIXIIIl 

Although we may observe the common tendency that the three form factors Fv{q'^), Friq'^) and FA{q^) become 
flatter as a function of q^ in the smallest volume (L — 12), we do not make a definite conclusion through this study. 
Rather we can say that the finite volume effects on the values of any form factor at finite q^ are less appreciable than 
our expectation raised by the fact that the axial- vector coupling significantly suffers from the finite volume effect. 
Indeed, we could attempt the direct comparison between results for two lattice sizes only at q^ « 0.44 GeV^, which 
is a relatively high value. Therefore, we deduce that the finite volume effects are not so serious in the high q^ region. 
Finally, it is worth mentioning that we confirm that our observed forms for the five form factors well describe the 
g^-dependence of those form factors even in the relatively high q^ region, up to at least « 1.0 GeV^, apart from 
consideration of finite volume effects. 



VI. COMPARISON WITH PREVIOUS RESULTS 



As we discussed in Sec. IIVI the larger spatial volume enables us to perform the shorter ^^-extrapolation to extract 
fundamental information on the nuclcon structure, e.g. the magnetic moment, the charge radius and the induced 
pseudo-scalar coupling from respective form factors without large systematic uncertainties. However, some of previous 
studies are performed on relatively small volumes, where the longer ^^-extrapolation is inevitable. In this context, 
our lattice setup is superior to previous studies. (See Table IXTVl for a summary of previous calculations 0, 0, SI3-) 

A large volume simulation, which is comparable to our lattice volume, (3.6 fm)'^, has been done by the LHPC 
Collaboration with the mixed action simulation using DWF valence quarks on the asqtad-improved gauge configura- 
tions with fourth-rooted staggered sea quarks. The quenching effect in our simulations could be observed through a 
comparison with the mixed action results. However, there is no detailed analysis of the usual form factors in Ref [H]- 
Instead, we can access their raw data of four nucleon form factors from their tables. We simply compare our measured 
form factors at the lightest quark mass (r7i7r=0.39 GeV) with their results with the lightest pion mass of 0.35 GeV 
and the largest volume of (3.5 fm)'^ in Fig. 1331 Surprisingly, all figures for four form factors show good consistency 
between our quenched results and the LHPC mixed action results within statistical errors, at least, in the low q^ 
region. This indicates that (un)quenching effects on these form factors are still small for ^ 0.35 GeV. However, 
this conclusion is rather premature since the mixed action simulation is not a fully dynamical simulation, rather a 
partially quenched simulation 51|]. We must wait for fully dynamical DWF simulation to make a firm conclusion. 
The RBC and UKQCD Collaborations have begun 2+1 flavor DWF calculations with large physical volume [s^. [s^. 
We will do such a comparison in a future publication. 

Finally, it is worth comparing our DWF results with the results obtained using Wilson fermions. In Table IXVl our 
results of the axial-vector coupling gA and the rms radii of the iso-vector Dirac and Pauli form factors, the iso-vector 
nucleon magnetic moment /Xp — /j,„, the axial dipole mass, the induced pseudo-scalar coupling gp and the pion-nucleon 
coupling g-n-NN in the chiral limit are compared with previous quenched Wilson results from Refs. Q and 0- Their 
lightest pion mass of 0.41 GeV and physical volume of (2.9 fm)"^ are relatively similar to our lattice set up. One finds 
that the rms radii and the axial dipole mass are quite consistent with each other, while the better agreement with the 
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experiment for the axial-vector coupling and the iso-vector nucleon magnetic moment appear in our quenched DWF 
results. Note that although the induced pseudo-scalar form factor Fp{q^) was calculated in Rcf. [7], the value of the 
induced pseudo-scalar coupling gp was not evaluated there. 

Let us estimate gp from their fit parameters given in Table III of Ref. 0^ as follows. First we read off app^ from their 
parameters obtained by the monopole fit co/(l-|-(?^/A^) to the ratio of 2M]s[Fp{q^)/ FAiq^)- The authors reported that 
the monopole mass A is larger than their measured pion mass, while the value of cq is smaller than 4M^ /'m-'i evaluated 
by the measured pion and nucleon masses. Therefore, their app^ may have appreciable g^-dependence, which is given 
by (1 + (7^/m^)/(l -I- q^ /h?). On the other hand, the value of app^ at = can be given by ComJ/(4M|f), which 
is ranged from 0.68 to 0.62 when the pion mass vary from 0.56 GeV to 0.41 GeV. There is the descending tendency 
with the decrease of the pion mass. A simple linear extrapolation with respect to the pion mass squared leads to a 
value in the chiral limit less than 0.6, that is significantly deviated from unity. Thus, one may easily deduce that 
their gp should be much smaller than our observed gp since their corresponding a-^-p^ and measured -Fa(O), both of 
which are main ingredients in Eq. (|48p. are about 40% and 20% smaller than our DWF values respectively. Indeed, 
our quoted of the quenched Wilson results in Table IXVl which is determined with the value of a^^^ evaluated at 

= 0.88m^, is almost a half of our quenched DWF value. This discrepancy is attributed to the fact that the Wilson 
fermions do not yield the correct pion-pole structure of Fp{q^) 0], while the PPD form provides a good description 
of the g^-dependence of Fp{q'^) at low q^ in our DWF calculation. 

VII. CONCLUSION 

In this paper, we have studied the weak nucleon form factors at low q^ in quenched lattice QCD. We have used 
domain wall fermions in a very large physical volume (3.6 fm)'^. There are two reasons for requiring such large 
volume. As shown in the early calculation of the axial- vector coupling gA p^ - the nucleon matrix element may suffer 
significantly large finite volume effects. However, we really did not know whether the spatial volume (2.4 fm)'^, that 
was utilized in Ref. [iB'l, was large enough for the nucleon. Secondly, the large spatial extent provides the capability 
to access lower non-zero momentum transfer. For the spatial volume (3.6 fm)'^, the smallest value of non-zero q^ is 
about 0.1 GeV^ 

We first demonstrated that the finite volume effect on the axial-vector coupling gA is well described by the power- 
law behavior, while the vector coupling gv has no appreciable finite volume effect. However, it is found that a serious 
finite volume effect on the axial-vector coupling gA is not seen in the range of the spatial lattice size from 2.4 fm to 
3.6 fm. Finally, we obtain the ratio of the axial-vector to the vector couphng [gA/gvY^^ = 1-219 ± 0.038(stat) ± 
0.024(norm) ± 0.002(vol) at the physical point from the largest volume (3.6 fm)'^, which agrees with our early estimate 
from the volume (2.4 fm)'^ and underestimates the experimental value of 1.2695(29) by less than 5%. 

Using the largest volume (3.6 fm)'^, we studied four of the weak nucleon form factors and also the pseu do-scalar form 
factor. The g^-dependences of all measured form factors at low q^ are discussed with great interest [64I. It is observed 
that the vector (i^y), induced tensor {Ft) and the axial- vector {Fa) form factors are well described by the dipole form 
as Fi{q^) = F.,{Q)I{1 + q^ /M^) {i = V,T,A) at low (g^ < 0.44 GeV^). Each measured dipole mass overestimates 
the corresponding experimental value by about 20 %. This fact indicates that corresponding rms radii are somewhat 
smaller than the actual nucleon size. However, interestingly, the ratios of dipole masses, Ma/Mv = 1.285 ±0.073 and 
Mt/Mv = 0.869 ± 0.057 are fairly consistent with respective experimental ones. We also calculated the difference of 
the proton and neutron magnetic moments, fJ,p — — 4.13 ± 0.23, from the value of FT{q^) at q^ = 0. Our obtained 
value is about 14% smaller than experimental value of 4.70589. 

We have presented a detailed study of the induced pseudo-scalar form factor Fp(g^), which is less well-known 
experimentally. It is observed that the g^-dependence of Fp{q^) exhibits the strong quark-mass dependence in the 
low q^ region. This is associated with the pion-pole contribution. Indeed, we confirm that the measured value of Fp{q^) 
is well described by the pion-pole dominance (PPD) form Fp™{q^) = 2MiqFA{q^) / {q^ + ?7iJ) with a multiphcative 
factor appD(< 1)- With the help of such the PPD-like form, we can evaluate the induced pseudo-scalar coupling as 
(ffp)™" = 8.15 ± 0.54(stat) ± 0.16(norm). This value is to be compared with the most recent experimental value of 
7.3 ± 1.1 from the MuCap experiment. Furthermore, we evaluated the pion-nucleon coupling g-„NN from the residue 
of the pion pole in the induced pseudo-scalar form factor and found gT^NN = 10.4 ± l.O(stat). 

We have also studied the axial Ward-Takahashi identity in terms of the nucleon matrix elements, which may be 
referred to as the generalized Goldberger-Treiman relation. For this purpose, we also calculated the pseudo-scalar 
matrix element, which is described by the single form factor called as the pseudo-scalar form factor Gp{q^). We 
have found that the measured g^-dependence of Gp{q^) is quite consistent with an expected behavior associated 
with FA{q'^) and Fp{q^) in consequence of the axial Ward-Takahashi identity, or the generalized Goldberger-Treiman 
relation. This fact ensures that the PPD form provides a good description of the g^-dependence of Gp{q'^) as well. 

In the case of large but finite fifth dimension, the axial Ward-Takahashi identity for domain wall fermions can be 
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slightly modified by introducing an additive shift of the quark mass due to the presence of the midpoint contribution 
to the divergence of the axial- vector current. Such an additive constant shift known as the residual quark mass 
is usually measured through mesonic two-point correlation functions. In an earlier calculation, the residual quark 
mass is observed much smaller than the lightest quark mass utilized here [20|. However, our observed quark mass 
shift rrishift, which is required for satisfaction of the generalized Goldberger-Treiman relation^_is close to 50% of our 
lightest quark mass and an order of magnitude larger than the residual mass quoted in Ref. [20]. This issue may be 
connected with the correct chiral behavior of the nucleon matrix elements. In the quenched approximation, there may 
be unknown quenching effects in nucleon matrix elements in the vicinity of the chiral limit. In this context, the above 
issue is beyond the scope of this quenched study. Rather, dynamical simulations are much preferable to investigate it 
further. The RBC and UKQCD Collaborations have begun Nf = 2-1-1 flavor domain wall fermion calculations with 
large physical volume V (2.7 fm)'^ and the lightest ud quark mass down to 1/7 the strange quark mass (m^ w 330 
MeV) [53| . We plan to develop the present calculation for a precise determination of nucleon form factors and also to 
address all of unsolved issues described in this paper. Such planning is now underway [ssj . 

Note added: After the completion of this work, we became aware of a paper [3] where the nucleon axial-vector 
form factor FAiq^) and induced pseudo-scalar form factor Fp[q^) are calculated in quenched and unquenched lattice 
QCD using Wilson fermions. 
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Appendix A: Various rms radii in the vector channel 

In Table [H the electric charge and magnetization radii for the proton and neutron are summarized. Using these 
experimental values, the iso- vector electric charge and iso- vector magnetization radii can be evaluated by the following 
relations 0, [1] 



where G^(m)(9^) 
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0.939(5) fm and 

\/((^m)^) = 0.862(14) fm. Similarly, the rms radii for the iso- vector Dirac form factor F^^q^) — Ffi^q^) — FlHq^) and 
the iso- vector Pauli form factor (9^) = -^2(1^) ~ P'li'f) be given through the following relations 0, 



E(M)^'f) ~ GE(Ai)iQ^) and ^ly = f^p - ^in- Then one obtains 
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which yield vaMP) = 0.797(4) fm and V(MF> = 0.879(18) fm. 



Appendix B: Generalized Goldberger-Treiman relation and pion-pole dominance 

The generalized Goldberger-Treiman relation is derived from the nucleon matrix elements of the currents on both 
sides of the axial Ward-Takahashi identity [l^]; daA'^{x) — 2rhP°-{x) where the exact iso-spin symmetry is considered 
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as m = m„ = m^. The nucleon matrix element of the divergence of the axial- vector current is represented in the 
following form: 

- [2MNFA{q^) - q^Fp{q^)]u,ip'h5t''uAp). (71) 

Here, it is worth mentioning that we have used the Dirac equation for the nucleon, t2„ (p)(if)+MN) — (ip+MN)uj^ (p) = 
to get from the first line to the second line. Then one easily finds that the g^-dependences of three form factors are 
constrained by the following relation 

2MjvF^(g') - g'Fp(g2) + 2mGp{q''), (72) 

which is a consequence of the axial Ward-Takahashi identity. This expression may be referred to as the generalized 
Goldberger-Treiman relation [2^1 • 

Here we discuss the case where the limits m and ^ are taken on Eg. ((72)) . Of course, the left-hand side 
(l.h.s.) of Eq. ([72|) yields a non-zero value in the double limit. First, we consider the case where the chiral limit is first 
taken before the limit oi q^ ^ 0. 



lim lim 2MNFAiq^) = lim lim Fp{q^) , (73) 

which requires the massless pion pole in Fp{q^) in the chiral limit [s^ as limm_^o Fp{q^) oc for non- vanishing of 
the l.h.s. of Eq. ([73]) . Secondly, the chiral limit is taken after the limit of q^ 0: 



lim lim 2MNFA(q^) ^ lim 2m lim Gp(q^)] , (74) 

which requires the 1/rh singularity in Gp{q^) at = as limq2_^Q Gp{q^) (x i ~ for non-vanishing of the l.h.s. 

of Eq. ll74l) . As a result, Fp{q^) and Gp{q^) must have the pion-pole structure which should become dominant at 
low q^ [33|. Therefore, one can deduce that Fp{q^) and Gp{q^) are described by the following forms, at least, in the 
vicinity of the pole position q^ = —rn^ [ssl. l49j. 

PPD 2 2MNFA{q^) 

Fp (<? ) = 5- (75) 

q'^ + 

2mGPPD(g2) = 2A/ArF^(g^) ^ ^ , (76) 

q'^ + 

which we call the pion-pole dominance (PPD) forms. Consequently, we realize that the ratio of G^^{q^) and 
Fp^^iq^) gives the low-energy constant Bq as 

where mi = 2rhBo. 
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In this paper, we restrict ourselves to considering the iso-spin symmetric case as m„ = md, where the second-class form 
factors do not appear in consequence of G-parity invariance [25[. 

The sign of all form factors is chosen to be positive. Remark that our 75 definition, 75 = 70,7^727* = —75^, has the opposite 
sign relative to that in the Minkowski convention (7*^ — i-y and 70^ = 74) adopted in the particle data group [22] . 
One might worry about the large scaling violation, which is observed for the axial-vector coupling qa in a quenched 
calculation with overlap fermions at the same lattice spacing [s^]. However, the previous quenched DWF studies reported 
that there is no appreciate scaling violation in the kaon B-parameter Bk [39| and proton decay matrix elements ]40| ] at 
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/3 = 0.87 (a ~ 0.15 fm) and 1.04 (a ~ 0.10 fm). We may deduce that no large scaling violation is ensured for other matrix 
elements as well in our DWF calculations. 
[57] We treat our data sets as 70 independent measurements after taking average of the multiple source results on each 
configuration. 

[58] It is worth mentioning that such consistency is not observed at L = 16, since the lowest point at L = 16 may suffer 

from the finite volume effect as we will describe in Sec. |V] 
[59] Here we use a traditional convention as = /tt/x/S ~ 93MeV, while is quoted in Ref . [20l| . 

[60] Strictly speaking, the relativistic dispersion relation (|18|) . which ensures that [if>+ Mm)Uj^(p) = 0, should be well satisfied 
in simulations. 

[61] In Ref. Q, the similar deviation from unity as a^wT > 1 '^^'^ read off from their Fig. 19 where an inverse of Q^wt 
given in the limit of — > 0. 

[62] Detailed knowledge about the g^-dependence of the weak form factors in neutron beta decay is required for our lattice 
study of flavor SC/(3) breaking effects in hyperon beta decays [50| . 
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TABLE I: Experimental values of magnetic moments, electric charge and magnetization radii of the proton and neutron. 



Observable Experimental value Reference 

lip +2.792847351(28) \22} 

fj,n -1.91304273(45) [22] 

{{r^ff^^ 0.8750(68) fm \22} 

{{rlf) -0.1161(22) fm^ [22] 

{{rljf)^^^ 0.855(35) fm [26] 

{[rl,f)^'^ 0.873(11) fm [27] 



TABLE IL Summary of available experimental data for the induce pseudo-scalar form factor Fp{q^). The smallest point is 
given by the MuCap experiment, while other three points are obtained from a single experiment of pion electroproduction 
at threshold. 



q' (GeV^) 


Fp(g2) (MeV"') 


Experiment (reference) 


0.0098 


0.069 ± 0.010 


ordinary muon capture [iL] 


0.073 


0.0229 ± 0.0028 


pion electroproduction [12. 34] 


0.139 


0.0140 ± 0.0022 


pion electroproduction [12. 341 


0.179 


0.00932 ± 0.00248 


pion electroproduction [12, 34] 



TABLE IIL Simulation parameters for each volume studied in this work. 



Gauge action (/3) 




X T 




Ms 


Quark mass values (m/) 


spatial size L [fm] 


Statistics 


# of sources 


DBW2 (0.87) 


24^ 


X 32 


16 


1.8 


0.02, 0.04, 0.06, 0.08 


3.6 


70 


3 




16« 


X 32 


16 


1.8 


0.04, 0.06, 0.08 


2.4 


377 


1 




12^ 


X 32 


16 


1.8 


0.04, 0.06, 0.08 


1.8 


800 


1 



TABLE IV: The residual mass rrires, inverse lattice spacing (flp ^, set by the p meson mass), the renormalization factor of the 
axial-vector current (Za), and the pion decay constant (^V). Those values are taken from Ref. [l^j, where simulations are 
performed on a 16^ X 32 volume. 



Gauge action (/3) 


Ms 


Ls 




ap' [GeV] 


ZA(mf = — rriros) 


[MeV] 


DBW2 (0.87) 


1.8 


16 


5.69(26) X 10"* 


1.31(4) 


0.77759(45) 


91.2(5.2) 
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TABLE V: Fitted masses of the pseudo-scalar meson state and fitted energies of the nucleon state with the five lowest momenta 
for each volume. All tabulated values are given in lattice units. Results for the nucleon energies with nonzero momenta are 
averaged over all possible permutations of the lattice momentum p = {na:,nx,nz) in units of 2n/L, including both positive and 
negative directions. 









En{p) 










size L 


ruf 




(0,0,0) 


(1,0,0) 


(1,1,0) 


(1,1,1) 


(2,0,0) 


24 


0.02 


0.3003 (10) 


0.8651 (70) 


0.9038 (79) 


0.9432 (94) 


0.9824 (117) 


1.0157 (130) 




0.04 


0.4143 (11) 


0.9801 (52) 


1.0152 (58) 


1.0496 (66) 


1.0829 (77) 


1.1161 (83) 




0.06 


0.5040 (11) 


1.0783 (47) 


1.1098 (50) 


1.1408 (55) 


1.1710 (60) 


1.2006 (66) 




0.08 


0.5829 (11) 


1.1690 (43) 


1.1976 (45) 


1.2257 (49) 


1.2531 (53) 


1.2793 (57) 


16 


0.04 


0.4148 (9) 


0.9869 (50) 


1.0595 (58) 


1.1251 (74) 


1.1822 (106) 


1.2436 (150) 




0.06 


0.5050 (8) 


1.0821 (42) 


1.1483 (48) 


1.2095 (61) 


1.2632 (84) 


1.3189 (111) 




0.08 


0.5837 (8) 


1.1703 (37) 


1.2313 (43) 


1.2886 (54) 


1.3391 (71) 


1.3890 (92) 


12 


0.04 


0.4150 (10) 


0.9795 (75) 


1.097 (10) 


1.190 (17) 


1.332 (68) 


1.236 (79) 




0.06 


0.5046 (9) 


1.0729 (55) 


1.1844 (70) 


1.280 (12) 


1.419 (41) 


1.466 (22) 




0.08 


0.5832 (8) 


1.1703 (37) 


1.2313 (43) 


1.289 (54) 


1.339 (71) 


1.389 (92) 



TABLE VI: Results for the vector coupling g^"'''", the axial-vector coupling gJ|-*"<^« and their ratio pj"^ = (sv/sa)'''"''"'- Gauss- 
smeared-to-gauss-smeared quark propagators are used in the present study, while box-to-local quark propagators were used in 
the previous calculation 15]. L — 16 results in the present calculations agree well with the previous L = 16 results. 



Smearing type (Ref ) 


X iVt 




(<?v)latticc 


(ffA)latticc 


(S'A)rcn 


Gauss-Gauss (this work) 


24^ X 32 


0.02 


1.2435(60) 


1.509(52) 


1.212(42) 






0.04 


1.2386(18) 


1.537(27) 


1.240(22) 






0.06 


1.2306(12) 


1.533(18) 


1.245(15) 






0.08 


1.2206(10) 


1.529(14) 


1.252(11) 




16^ X 32 


0.04 


1.2429(24) 


1.494(29) 


1.202(24) 






0.06 


1.2332(12) 


1.497(17) 


1.214(14) 






0.08 


1.2224(9) 


1.499(13) 


1.226(11) 




12^ X 32 


0.04 


1.2465(34) 


1.441(51) 


1.156(41) 






0.06 


1.2328(16) 


1.462(27) 


1.185(22) 






0.08 


1.2213(11) 


1.474(17) 


1.206(14) 


Box-Local fl5] 


16^ X 32 


0.02 


1.2440(30) 


1.531(60) 


1.229(49) 






0.04 


1.2323(14) 


1.523(24) 


1.230(20) 






0.06 


1.2220(10) 


1.510(15) 


1.230(12) 






0.08 


1.2106(8) 


1.505(11) 


1.236(9) 




8^ X 24 


0.04 


1.223(10) 


1.303(146) 


1.059(120) 






0.06 


1.214(5) 


1.342(74) 


1.099(62) 






0.08 


1.203(4) 


1.373(46) 


1.136(39) 
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TABLE VII: Fitted results of Fi(g^) and -^2(9^) with the dipole form [Eqs. ((37|l and ((39])] and their extrapolated values to the 
chiral limit and the physical point. 





Fi(g2) = F^(g2) 


{rl)i (fm) 


F2(g2) = 2A/jvi?T(g') 






TUf 


Mv (GeV) 


FriO)=F2iO)/Fv{0) 


Mt (GeV) 


(4)' (fm) 


0.08 


1.330(17) 


0.514(7) 


3.76(7) 


1.155(18) 


0.592(9) 


0.06 


1.286(21) 


0.532(9) 


3.65(9) 


1.112(22) 


0.615(12) 


0.04 


1.239(32) 


0.552(14) 


3.48(13) 


1.072(32) 


0.637(19) 


0.02 


1.177(65) 


0.581(32) 


3.08(26) 


1.062(78) 


0.644(47) 


phys. point 


1.148(55) 


0.589(25) 


3.13(23) 


1.002(58) 


0.676(34) 


chiral limit 


1.142(57) 


0.592(25) 


3.11(24) 


0.997(60) 


0.679(35) 


Empirical values 


0.857(8) 


0.797(4) 


3.70589 


0.778(23) 


0.879(18) 



TABLE VIII: Fitted results of Fa(i3^) with the dipole form [Eq.(l37ll] and Fp(g2) with the PPD-like form [Eq.(|47l)], and their 
extrapolated values to the chiral limit and the physical point. 











Fp{q') 




nif 


(3,4)^=" = Fa(0)/F^(0) 


A^A (GeV) 


{rl}i (fm) 


(ffp)^^" 




0.08 


1.252(11) 


1.572(26) 


0.435(7) 


11.08(28) 


0.815(18) 


0.06 


1.245(15) 


1.541(32) 


0.444(9) 


10.37(28) 


0.831(18) 


0.04 


1.240(22) 


1.523(47) 


0.449(14) 


9.66(34) 


0.853(20) 


0.02 


1.212(42) 


1.618(126) 


0.422(33) 


8.73(51) 


0.884(31) 


phys. point 


1.219(38) 


1.502(83) 


0.453(24) 


8.15(54) 


0.897(32) 


chiral limit 


1.218(40) 


1.500(85) 


0.454(25) 


8.04(55) 


0.900(33) 


Experimental values 


1.2695(29) 


1.026(21) 


0.666(14) 


7.3 (1.1) 


1 (Theor.) 



TABLE IX: Ratio of the dipole masses and their extrapolated values to the chiral limit and the physical point. 



ruf 


Mt I My 


MaIMv 


0.08 


0.869(14) 


1.182(17) 


0.06 


0.865(19) 


1.198(24) 


0.04 


0.865(32) 


1.229(42) 


0.02 


0.902(84) 


1.374(131) 


phys. point 


0.869(57) 


1.285(73) 


chiral limit 


0.869(58) 


1.289(75) 


Empirical values 


0.908(28) 


1.197(27) 



TABLE X: Results of the ratio a^wTi ^^e quark mass (jtIawt) defined in the generalized Goldberger-Treiman relation [Ea. (|62|l ] 
and the mass shift mshift. 



m/ 




"^AWT 


mshift = "^AWT ^ 


0.08 


1.12(1) 


0.0896(7) 


0.0096(7) 


0.06 


1.15(1) 


0.0688(7) 


0.0088(7) 


0.04 


1.21(2) 


0.0484(8) 


0.0084(8) 


0.02 


1.41(6) 


0.0281(12) 


0.0081(12) 





N/A 


0.0073(12) 


0.0073(12) 
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TABLE XI: Results of five (dimensionless) form factors computed on a 24^ x 32 volume. All form factors are renormalized 
except for the pseudo-scalar form factor Gp{<f). For notational simplicity, we use the momentum q in units of 2-k/L. 



ruf 


q 


(GeV^) 




2A/ivJ'^.™(g') 


Fr{q^) 


2AfjvFF"(g') 


G),^'%q^) 


0.02 


(0,0,0) 


0.000 


1.0000(48) 


N/A 


1.212(42) 


N/A 


N/A 




(1,0,0) 


0.113 


0.852(15) 


2.53(21) 


1.132(41) 


19.79(1.55) 


21.75(1.11) 




(1,1,0) 


0.222 


0.739(23) 


2.16(17) 


1.040(40) 


12.00(99) 


15.26(79) 




(1,1,1) 


0.326 


0.633(29) 


1.82(17) 


0.892(40) 


8.36(74) 


11.76(78) 




(2,0,0) 


0.427 


0.617(42) 


1.61(19) 


0.921(62) 


7.40(79) 


8.88(1.01) 


0.04 


(0,0,0) 


0.000 


1.0000(15) 


N/A 


1.240(22) 


N/A 


N/A 




(1,0,0) 


0.114 


0.864(6) 


2.88(11) 


1.131(21) 


15.61(81) 


17.43(43) 




(1,1,0) 


0.224 


0.757(10) 


2.42(10) 


1.031(21) 


10.97(50) 


13.12(34) 




(1,1,1) 


0.331 


0.668(13) 


2.06(9) 


0.927(19) 


8.20(43) 


10.45(33) 




(2,0,0) 


0.434 


0.622(18) 


1.83(10) 


0.897(28) 


6.91(47) 


8.44(35) 


0.06 


(0,0,0) 


0.000 


1.0000(10) 


N/A 


1.245(15) 


N/A 


N/A 




(1,0,0) 


0.114 


0.872(4) 


3.05(8) 


1.133(14) 


13.46(59) 


14.66(25) 




(1,1,0) 


0.225 


0.771(6) 


2.60(7) 


1.036(14) 


10.27(39) 


11.60(20) 




(1,1,1) 


0.333 


0.687(8) 


2.24(7) 


0.946(13) 


8.01(33) 


9.53(20) 




(2,0,0) 


0.439 


0.633(11) 


1.98(7) 


0.895(18) 


6.72(39) 


7.87(22) 


0.08 


(0,0,0) 


0.000 


1.0000(8) 


N/A 


1.252(11) 


N/A 


N/A 




(1,0,0) 


0.114 


0.879(3) 


3.18(6) 


1.143(10) 


12.29(53) 


12.88(17) 




(1,1,0) 


0.226 


0.782(4) 


2.74(6) 


1.048(10 


9.80(34) 


10.51(14) 




(1,1,1) 


0.335 


0.702(6) 


2.39(5) 


0.963(10) 


7.88(29) 


8.82(14) 




(2,0,0) 


0.442 


0.645(8) 


2.11(6) 


0.905(14) 


6.67(37) 


7.41(17) 



TABLE XIL The same as Table |Xll but for a 16^ x 32 volume. 



rUf 


q 


q^ (GeV^) 


F{F"(g2) 


2MivF^<="(g2) 




2MjvFi,™(g2) 


G^="^'=(g^) 


0.04 


(0,0,0) 


0.000 


1.0000(19) 


N/A 


1.202(24) 


N/A 


N/A 




(1,0,0) 


0.251 


0.718(15) 


2.14(10) 


0.939(23) 


7.87(64) 


12.54(45) 




(1,1,0) 


0.485 


0.565(19) 


1.53(8) 


0.751(30) 


6.00(40) 


7.45(36) 




(1,1,1) 


0.706 


0.430(25) 


1.14(8) 


0.645(35) 


3.78(34) 


5.04(35) 




(2,0,0) 


0.915 


0.417(49) 


0.97(14) 


0.558(64) 


2.29(43) 


4.53(59) 


0.06 


(0,0,0) 


0.000 


1.0000(10) 


N/A 


1.214(14) 


N/A 


N/A 




(1,0,0) 


0.253 


0.744(9) 


2.37(7) 


0.963(15) 


8.00(52) 


10.93(25) 




(1,1,0) 


0.491 


0.593(13) 


1.80(5) 


0.788(20) 


6.04(27) 


7.09(21) 




(1,1,1) 


0.717 


0.471(17) 


1.33(6) 


0.688(23) 


4.17(25) 


4.93(22) 




(2,0,0) 


0.933 


0.430(29) 


1.06(9) 


0.580(37) 


2.48(32) 


4.15(34) 


0.08 


(0,0,0) 


0.000 


1.0000(7) 


N/A 


1.226(11) 


N/A 


N/A 




(1,0,0) 


0.254 


0.760(7) 


2.53(5) 


0.988(11) 


7.99(50) 


9.92(18) 




(1,1,0) 


0.495 


0.611(10) 


1.93(4) 


0.816(15) 


6.15(22) 


6.76(15) 




(1,1,1) 


0.725 


0.496(13) 


1.47(5) 


0.714(17) 


4.42(22) 


4.83(16) 




(2,0,0) 


0.946 


0.440(20) 


1.14(7) 


0.599(27) 


2.64(30) 


3.98(25) 
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TABLE XIII: The same as Table IXB but for a 12^ x 32 volume. 



771/ 


q 


(GeV^) 




2A/ivF^<="(g^) 




2MjvFF"(g') 




0.04 


(0,0,0) 


0.000 


1.0000(27) 


N/A 


1.156(41) 


N/A 


N/A 




(1,0,0) 


0.434 


0.681(34) 


2.02(14) 


0.883(43) 


7.47(74) 


7.81(61) 




(1,1,0) 


0.822 


0.466(53) 


1.32(16) 


0.693(81) 


3.44(56) 


4.79(66) 


0.06 


(0,0,0) 


0.000 


1.0000(13) 


N/A 


1.185(22) 


N/A 


N/A 




(1,0,0) 


0.439 


0.672(17) 


2.03(8) 


0.885(23) 


7.00(54) 


7.60(33) 




(1,1,0) 


0.837 


0.468(25) 


1.33(8) 


0.666(38) 


3.59(34) 


4.41(30) 


0.08 


(0,0,0) 


0.000 


1.0000(9) 


N/A 


1.206(14) 


N/A 


N/A 




(1,0,0) 


0.442 


0.671(11) 


2.09(6) 


0.891(16) 


6.74(48) 


7.23(23) 




(1,1,0) 


0.848 


0.477(15) 


1.38(5) 


0.669(25) 


3.74(25) 


4.20(19) 



TABLE XIV: Previous lattice calculations for nucleon form factors. Note that our lowest non-zero (gmin) is smaller than 
previous calculations. This is an essential point in determinations of the induced pseudo-scalar coupling {gp) from Fp{q^) 
and the nucleon magnetic moments {(j-n) from Friq'^) in order to reduce the systematic uncertainties stemming from long 
^^-extrapolation. 



Group (reference) 


Type 


Fcrmion (valence) 


a [fm] 


spatial size [fm] 


[GeV] 


qLn [GoV^] 


QCDSF [2] " 


Quench 


Clover 


0.11 


1.8 


0.54, 0.64, 0.74, 0.91, 0.98 


0.47 








0.08 


1.9 


0.61, 0.76, 0.89, 1.03 


0.40 








0.06 


1.9 


0.63, 0.79, 0.79, 0.93, 1.05 


0.39 


Cyprus-MIT [4, 7] 


Quench 


Wilson 


0.09 


2.9 


0.41, 0.49, 0.56 


0.17 




FuU {Nf = 2) 


Wilson 


0.08 


1.9 


0.38, 0.51,0.69 


0.42 


LHPC [5] ' 


Mixed {Nf =2 + 1) 


DWF 


0.124 


3.5 


0.35 


0.11 








0.124 


2.5 


0.36, 0.50, 0.60, 0.68, 0.76 


0.18-0.20 


This work 


Quench 


DWF 


0.15 


3.6 


0.39, 0.54, 0.66, 0.76 


0.11 








0.15 


2.4 


0.54, 0.66, 0.76 


0.25 








0.15 


1.8 


0.54, 0.66, 0.76 


0.43 



"Only electric-magnetic form factors {Fy, Ft) are studied. 

'Raw data of Fy, F^, i^A and Fp are available in tables, while detail analysis is not found. 

'^DWF valence quarks on the asqtad-improved gauge configurations with fourth-rooted staggered sea quarks. 



TABLE XV: Comparisons with previous quenched Wilson results obtained from Refs. [j and Q. The extrapolated values in 
the chiral limit are evaluated by a simple linear quark mass dependence. 



Reference 


f-l'P f^rt 


{rl}i [fm] 


(4)^ [H 


(Sa)-^^" 


Ma [GeV] 


(ffp)'^" 




Refs. [4] and [7] 


3.73(13) 


0.585(13) 


0.72(2) 


1.065(24) ' 


1.500(59) * 


3.54(61) " 


11.8(0.3) 


This work 


4.11(24) 


0.592(25) 


0.679(35) 


1.218(40) 


1.500(85) 


8.04(55) 


10.4(1.0) 



"g-wNN was evaluated at = in Ref. [7j, while our Qt^nn is defined at the pion pole, = —m^ with the physical pion mass. 
^ The chiral extrapolated values are not available in Ref [7j. Therefore, we performed correlated fits with their measured data to evaluate 
them. 

" Although Fp(<j2) was calculated in Ref. the value of the induced pseudo-scalar coupling was not evaluated. The quoted value is 
estimated by us as described in the text. 
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FIG. 1: Comparison of measured and estimated energies of the nucleon for ruf = 0.04 as a function of absolute value of three- 
momentum jp|. Open circles, squares and diamonds, which corresponds to the measured values in lattice units, are obtained 
from L = 24, L = 16 and L — 12. The estimated energies are given by the relativistic dispersion formula -B(p) — -y^p^ + Mf, 
for continuum-like momenta — (dashed-dotted curve) and lattice momenta pi — sin[^ni] (dashed curve) with the rest 

mass Afjv measured ai L — 24. 
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FIG. 2: Relevant ratios of three- and two-point functions, (top) and Ay (bottom), for all possible three-momentum transfer 
q as a function of the current insertion time slice at m/ = 0.04. 
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FIG. 3: Relevant ratios of three- and two-point functions, Af;{qz — 0) (top), A^(g2 7^ 0) (middle) and Ay (bottom), for all 
possible three-momentum transfer q as a function of the current insertion time slice at m/ = 0.04. 
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FIG. 4: The physical ratio of couplings qa/qv as a function of the pion mass squared. Results on the largest volume (3.6 fm)^ 
(circles) exhibit milder quark mass dependence, while the smaller volume results (right-oriented triangles) show a slow downward 
tendency toward the chiral limit away from the experimental point (asterisk) . 
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FIG. 5: The vector coupling (gv')'^""^'' and axial-vector coupling (jta)'''**''^'' as functions of spatial lattice size for mf = 0.04 
(left figure) and m/ = 0.08 (right figure). Dashed curves are fits of the form ^^"'"'(■t') = gT^'^^oo) + bL'^. 



31 



0.9 



O 0.8 
> 

"~ 0.7 



1 ' ' ' 1 ' •- 

' \ 

\ 

\ 



0.5 - 



-1 1 1 1 r- 



m=0.02 
mj=0.04 



I 



^ ni=0.06 

A m=0.08 
— ■ experiment 



0.2 0.3 

q'[GeV'] 



FIG. 6: The renormalized Dirac form factor, F^'^{q^) = Fv{q^)/Fv(fi) as a function of four-momentum squared q 
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FIG. 7: Effective dipole-mass plot for the Dirac form factor Fv(q^) as a function of four-momentum squared q^ at m/ — 0.02. 
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FIG. 8: Chiral extrapolation of the Dirac dipole mass My- The extrapolated points in the chiral limit and at the physical 
point are represented by an open diamond and a filled diamond. The experimental value is marked with an asterisk. 
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FIG. 9: The renormalized Pauli form factor, i*2'^"(g^) — 2AIFT{q^)/Fv{0) as a function of four-momentum squared 
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FIG. 10: The extrapolation toward = for either F^^^Cq^) or GM{q'^) /Gsiq^) - 1 for nij = 0.02. Those quantities should 
intersect each other at — 0. For F2'^'^{q'^), the dipole form is applied, while a simple linear extrapolation with respect to q^ 
is used for Ghiiq^) / Gsiq^) — 1 thanks to its mild g^-dependence. The extrapolated values from both determinations agree well 
with each other, although those underestimate the experimental value. 
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FIG. 11: Chiral extrapolation of /ip — The square symbols have been moved slightly in the plus x-direction. 
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FIG. 12: Chiral extrapolation of the Pauli dipole mass Mt- Symbols are defined as in Fig. |8l 
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FIG. 13: The axial- vector form factor normalized by Fa{0) as a functiorr of four- momentum squared g^. 
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FIG. 14: Effective dipole-mass plot for the axial- vector form factor FA{q'^) as a function of four-momentum squared at 
m/ = 0.02. 
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FIG. 15: Chiral extrapolation of the axial dipole mass Ma- Symbols are defined as in Fig. |8] 




FIG. 16; Ratios, Ma/Mv and Mt /Mv as functions of squared pion mass. The extrapolated points in the chiral limit and at 
the physical point are represented by an open diamond and a filled diamond. The experimental values are marked with asterisk 
symbols. 
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FIG. 17: The renormalized and dimensionless induced pseudo-scalar form lactor 2Mjv-Fi?"(g^) as a function of four-momentum 
squared q^. 
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FIG. 18: The ratio of FF^C*?^) and as a function of four-momentum squared q 
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FIG. 20: The induced pseudo-scalar coupling {gpY'^'^ evaluated through two determinations, where measured values (squares) 
and the physical value (circles) are used for the nucleon mass in Ea.(|48p. 
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FIG. 21: The bare pseudo-scalar form factor Gp{q^) as a function of four-momentum squared 
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FIG. 22: The ratio Q^wt defined in Ea. (|59|l as a function of four-momentum squared 
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FIG. 23: The modified ratio rnf{aj^^^ ~ 1), which may be expressed as mAWT — nif, as a function of four-momentum squared 




FIG. 24: The mass shift mshift ~ "^awt ~ "^Z as a function of pion mass squared. For a comparison, the value of rrircs is 
included as a asterisk symbol. 
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FIG. 25: The ratio of Gp(Q^)/F|f"(g^) as a function of four-momentum squared q^. A slight q^-dependence remains against 
the naive expectation from the pion-pole dominance hypothesis. 
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FIG. 26: The ratio of Gp(g^)/FF"('?^)/AppD((7^) as a function of four- momentum squared . The residual g'^-dependence in 
the ratio of G p {q'^) / F^'^ (q'^) disappears by a multiplication of l/AppD(g^), which is responsible for the fact of appD 7^ 1. 
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FIG. 27: The ratio of Gp (0)/FF"(0)/Appd(0), m^/2m/, and m^/2(m/ + m.es) as functions of m'i. 
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FIG. 28: The normalized vector form factor Fv{q^) / Fv{0) obtained from simulations on lattices with three different spatial 
sizes. The left (right) panel is for ruf — 0.04 (m/ = 0.08). Filled circles, squares and diamonds are obtained from L — 24, 
L — 16 and L = 12. The curves represent the dipole form fits on results of the largest volume {L — 24). 
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FIG. 29: The renormalized and dimensionless induced-tensor form factor -F'2™(9'^) = 2MjvFT(g^)/-fV(0) obtained from simu- 
lations on lattices with three different spatial sizes. The left (right) panel is for m/ = 0.04 (m/ — 0.08). Symbols and solid 
curves are defined as in Figs [28] 
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FIG. 30: The normalized axial- vector form factor F_a((?^)/Fa(0) obtained from simulations on lattices with three different 
spatial sizes. The left (right) panel is for ruf — 0.04 (m/ = 0.08). Symbols and solid curves are defined as in Figs l28l 
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FIG. 31: The renormalized and dimensionless induced-pseudo-scalar form factor, 2A/jvJp°"(g^) obtained from simulations on 
lattices with three different spatial sizes. The left (right) panel is for nif — 0.04 (m/ — 0.08). Symbols are defined as in Figs l28l 
Solid curves are resulting fits with the form (|47|l on results of the largest volume {L = 24). 
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FIG. 32: The bare pseudo-scalar form factor Gp{<f') obtained from simulations on lattices with three different spatial sizes. 
The left (right) panel is for mf — 0.04 (m/ — 0.08). Symbols are defined as in Figs l28l Solid curves are resulting fits with the 
form (|63p on results of the largest volume (L = 24). 
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FIG. 33: Comparisons to results obtained from the LHPC mixed action results [^. Results for all four form factors at low q'^ 
are consistent with each other. This suggest that unquenching effects on these form factors are still small for ^ 0.35 GeV. 



